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ABSTBACT 

Two  niethods  for  evaluating  the  accuracy  of  hydrographic 
positioning  data  are  presented.   One  method  consists  of 
classifying  each  position  in  a  survey  based  on  the  radius  cf 
the  90  percent  confidence  circle.   The  second  method 
involves  classification  of  positions  based  on  the  parameters 
of  the  90  percent  confidence  ellipse.   Both  methods  are 
based  on  geometric  and  statistical  relationships  between 
intersecting  lines  of  position. 

Eange-range,  azimuth-azimuth,  and  range-azimuth  posi- 
tioning data  are  classified  using  both  criteria.   For 
noncritical  positions,  the  confidence  circle  method  is  found 
to  be  preferable  due  to  its  ease  of  interpretation.   Fox 
positions  of  significant  features,  such  as  underwater 
hazards,  the  confidence  ellipse  provides  a  more  useful 
representation  of  the  shape  and  orientation  of  the  true 
error  distribution. 

The  concept  of  presurvey  positioning  design  is  also 
presented.   With  the  aid  of  computer  graphic  displays,  the 
hydrographer  can  predict  the  accuracy  of  offshore  posi- 
tioning data  prior  to  data  acquisition.   By  analyzing  accu- 
racy lobes  generated  about  shore  stations,  a  survey  can  be 
designed  to  meet  given  specifications. 
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I.  INTRODUCTION 

2.   BACKGROUND 

A  hydrographic  record  can  he  viewed  as  the  resultant  cf 
two  independent  measurements  made  at  a  discrete  point  ever  a 
tody  of  water.   These  measurements  involve  the  determination 
of  a  vessel's  position  at  a  given  time  as  well  as  the  depth 
of  water  at  that  position.   Of  interest  to  the  hydrographer 
and  tc  the  user  of  hydrographic  data  is  the  accuracy  of  the 
position  determina tiers.   Fundamental  to  the  determination 
cf  positional  accuracy  is  the  identification  of  the  sources 
of  errors  in  positicn  measurements  and  the  ultimate  treat- 
ment cf  these  errors. 

A  hydrographic  position  can  be  determined  by  a  number  of 
methods  all  involving  geometric  relationships  between  known 
points  and  the  vessel's  unknown  location.   The  known  points 
may  be  fixed  stations  on  shore,  whose  coordinates  have  teen 
determined  by  geodetic  survey  methods,  or  they  may  be 
rapidly  iroving  satellites  whose  coordinates  in  time  and 
space  can  he  defined  very  precisely.   A  hydrographic  posi- 
tion is  established  by  the  intersection  of  two  or  more  lines 
of  position  (LOP's)  which  are  generated  by  the  geometric 
relationships  between  the  fixed  points  and  the  vessel's 
unknown  location.   The  resultant  accuracy  of  the  vessel's 
position  is  therefore,  in  part,  a  function  of  the  errors 
associated  with  the  intersecting  LOP's. 

Several  measures  cf  accuracy  can  be  used  to  evaluate  the 
guality  of  a  hydrographic  position.   Predictability,  or 
absolute  accuracy,  is  the  measure  of  accuracy  with  which  the 
positioning  system  can  define  the  location  of  the  same  point 
in  terms  of  geographic  coordinates.   Eepeatability ,  or 


relative  accuracy,  is  a  measure  with  which  a  positioning 
system  permits  a  user  to  return  to  a  specific  point  on  the 
earth's  surface  in  terms  of  the  LOP's  generated  ty  the 
system  [Bef.  1,  p.  1  4  ]-   With  the  elimination  of  all  system- 
atic cr  tias  errors,  the  terms  repeatability  and  predicta- 
bility become  identical.   Hydrcgraphic  surveyors  usually 
work  toward  this  condition,  although  it  is  not  always 
achievable. 

Heinzen  [Ref.  2]  and  Burt  [fief.  3]  have  presented 
several  techniques  for  quantifying  the  repeatable  accuracy 
for  offshore  positions.   These  techniques  have  roots  in  tie 
statistical  treatment  of  random  error.   Although  the  methods 
have  teen  well  documented,  no  single  criterion  to  classify 
the  accuracy  of  a  hydrographic  position  has  been  agreed  upon 
by  the  international  hydrographic  community. 

Preceding  the  development  of  automation  in  hydrographic 
data  acquisition  and  processing,  the  task  of  calculatirg  an 
accuracy  figure  to  attach  to  each  position  in  a  hydrographic 
survey  was  unthinkable.   To  ensure  overall  accuracy  in  a 
survey,  certain  generalizations  were  developed  to  act  as 
guidelines.   For  example,  the  U.S.  Coast  and  Geodetic  Survey 
Hydrographic  Manual  [Eef.  1,  p.  217]  states  the  following 
concerning  the  strength  of  a  three-point  fix: 

The  fix  is  strong  when  the  sum  of  the  two  angles  is 
equal  to  or  greater  than  180°  and  neither  angle  is  less 
than  3C°.   The  nearer  the  angles  equal  each  other  the 
stronger  will  be  tie  fix. 

Generalizations  of  this  type  provided  useful  qualitative 
guidance  for  assuring  a  degree  of  positional  accuracy  and 
many  are  still  in  existence  today. 

With  the  aid  of  computers,  the  hydrographer  now  has  the 
capacity  to  evaluate  the  accuracy  of  positioning  data  for  an 
entire  survey.   An  accuracy  figure  can  be  computed  for  each 
position  in  a  survey  and  stored  in  a  data  base  along  with 
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ether  survey  information.   This  figure  may  provide  useful 
information  for  users  of  the  data,  as  well  as  a  yardstick 
for  the  hydrographer  to  evaluate  the  quality  of  the  work. 
Furthermore,  a  presurvey  accuracy  analysis  enables  a  survey 
to  be  designed  to  meet  desired  specifications. 

E.   ACCOEACY  STANDARDS  FOB  HYDBOGEAPHIC  POSITIONING 

In  1S82,  the  International  Hydrographic  Organization 
(IHO)  published  new  recommendations  for  error  standards 
concerning  the  accuracy  of  hydrographic  positions.   Ihese 
standards  £Eef.  5]  are: 

The  position  of  soundings,  dangers  and  all  other  signif- 
icant features  should  be  determined  with  an  accuracy 
such  that  any  probable  error,  measured  relative  to  shore 
control,  shall  seldom  exceed  twice  the  minimum  plottahle 
error  at  the  scale  cf  the  survey  (normally  1.0  mm  on 
paper) .   It  is  most  desireable  that  whenever  positions 
are  determined  by  the  intersection  of  lines  of  position, 
three  such  lines  be  used.   The  angle  between  any  pair 
should  not  be  less  than  30°. 

Most  statisticians  define  the  term  "probable  error"  as 
that  errcr  occurring  at  the  50  percent  probability  level. 
However,  the  author  cf  the  IHO  standards,  Commodore  A.H. 
Cooper  RAN  (Ret.)  has  stated  that  the  term  "probable  error" 
was  intended  to  have  no  statistical  significance.   Munscn 
interpreted  the  words  "shall  seldom  exceed"  to  mean  10 
percent  of  the  time  £Bef.  6].   Using  this  interpretation, 
the  first  sentence  of  the  specification  might  be  written: 

The  position  of  soundings,  dangers  and  all  other  signif- 
icant features  should  be  determined  with  an  accuracy 
such  that  any  error  in  position  measured  relative  to 
shore  control  will  fall  within  a  circle  with  radius  of 
the  iririmum  plottable  error  at  the  scale  of  the  survey 
(normally  1.0  mm.  on  paper),  with  90  percent  confidence. 

The  specification  in  this  form  could  be  evaluated  quantita- 
tively.  The  criterion  for  defining  accuracy  in  terms  of  a 
fixed  probability  is  common  in  the  field  of  surveying.   For 
example,  the  standards  of  accuracy  developed  for  geodetic 
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control  surveys  have  their  origin  in  probability  thecry. 
Procedur€s  for  obtaining  first-crder  geodetic  positions 
require  sixteen  repeated  theodolite  observations  of  each 
direction.   Lower  order  positions  require  fewer  numbers  of 
observations.   Given  the  precision  of  one  observation  cf 
each  direction,  it  can  be  demonstrated  that  increasing  the 
number  of  observations  coincides  with  increasing  the  prob- 
ability cf  the  direction  falling  within  specified  limits. 
Begarding  accuracy  determinations,  there  are  several 
problems  unique  to  hydrographic  surveying.   Whereas  stan- 
dards for  other  types  of  surveys  rely  on  multiple  observa- 
tions of  the  same  quantity,  the  accuracy  of  a  hydrographic 
position  must  be  evaluated  in  terms  of  a  single  observation 
(which  may  be  the  intersection  of  two  or  more  LOP's). 
Diverse  methods  for  obtaining  a  hydrographic  position  exist 
and  these  methods  must  all  be  evaluated  using  the  same 
criterion.   Also,  there  is  a  broad  spectrum  of  equipment 
used  in  hydrographic  positioning  and  in  many  cases  the 
precision  of  this  eguipment  is  not  well  defined. 

C.   CEJECTIVES 

A  need  exists  to  give  quantitative  meaning  to  the  accu- 
racy specifications  set  forth  by  the  IHO.   One  of  the  objec- 
tives of  this  thesis  is  to  demonstrate  that  defining  the 
specifications  in  terns  of  the  fixed  90  percent  confidence 
level  is  a  valid  interpretation.   By  defining  what  the  spec- 
ifications imply,  procedures  can  be  developed  to  meet  the 
standards. 

A  second  objective  of  this  thesis  is  to  apply  tie  thecry 
of  errors,  associated  with  hydrcgraphic  positioning,  to  a 
data  set.   This  analysis  involves  classifying  positioning 
data  acquired  in  a  survey  based  on  the  radii  of  circles  of 
equivalent  probability.   It  will  be  demonstrated  that  this 
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method  of  classification  is  a  useful  index  for  quantifying 
the  accuracy  of  positions.   The  computed  radii  of  the  9C 
percent  confidence  circles  can  serve  as  an  accuracy  figure 
that  can  be  attached  to  each  position  in  a  survey  and  stored 
in  a  cata  base. 

The  third  objective  of  this  thesis  is  to  demonstrate 
that  a  presurvey  analysis  can  be  used  in  designing  posi- 
tional accuracy  to  meet  specifications.   The  existing 
general  guidelines  for  planning  can  be  better  defined.   For 
example,  in  planning  a  survey  hydrographers  usually  lay  out 
circles  hhich  delimit  the  30°  and  150°  boundaries  that 
define  the  minimum  and  maximum  allowable  intersectior  angles 
between  two  L0P*s.   As  a  means  to  meet  accuracy  require- 
ments, it  can  be  shown  that  these  limits  should  vary  based 
on  the  scale  of  the  survey  and  the  precision  of  the  posi- 
tionicg  eguipment. 
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II.  MATURE  OF  THE  PROBLEM 

The  development  cf  an  accuracy  figure  for  offshcre  posi- 
tions is  inherently  tied  to  the  geometry  of  the  positioning 
method  and  the  errors  which  are  associated  with  the  posi- 
tioning eguipment  that  is  used.   This  chapter  will  discuss 
the  geometric  and  statistical  elements  involved  in  deter- 
nining  an  offshore  position  and  presents  several  methods  for 
quantifying  repeatable  accuracy. 

A.   HYDECGRAPHIC  POSITIONING  GEOMETRIES 

An  offshore  fix  can  be  determined  by  the  intersection  cf 
two  or  more  LOP's.   These  LOP's  may  be  generated  by  elec- 
tronic or  visual  means.   Working  toward  the  development  of 
an  accuracy  index,  it  will  be  necessary  to  compute  the  angle 
of  intersection  of  tte  LOP's  associated  with  different  posi- 
tioning geometries.   The  following  sections  discuss  the 
geometry  cf  conventicnal  offshore  positioning  methods  and 
ways  to  compute  the  angles  of  intersection.   This  thesis 
will  not  address  the  geometry  involved  in  a  three-point 
sextant  fix. 

1 .   Eange-Ranqe 

Establishing  an  offshore  fix  by  range-range  geometry 
involves  measuring  distances  electronically  from  fixed  posi- 
tions on  shore  to  the  vessel's  unknown  location.   Ranges  can 
be  determined  by  measuring  the  elapsed  time  between  trans- 
mission and  receipt  of  a  radio  pulse  or  by  comparing  the 
phase  of  the  transmitted  wave  with  the  phase  of  the  received 
wave  [Ref.  2].   In  each  case,  transmitters  are  set  on 
stations  on  shore  whose  coordinates  are  determined  by 
precise  land  survey  irethods. 
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An  electronic  positioning  system  may  be  active  cr 
passive.   In  an  active  system,  a  transmitter  from  the  survey 
launch  keys  the  transmission  of  ranges  from  the  shore 
station.   In  turn,  the  signals  generated  from  the  shcre 
stations  (slaves)  are  then  received  by  the  launch.   An 
active  system  is  limited  to  a  finite  number  of  users, 
usually  not  more  than  about  four.   The  number  of  users  cf  a 
passive  system  is  unlimited  as  the  survey  launch  requires 
only  a  receiver  which  is  constantly  listening  for  signals 
which  are  being  transmitted  from  shore. 

Short-range,  cr  line-of-sight,  positioning  systems 
are  used  for  nearshore  hydrographic  surveys.   These  systems 
operate  in  the  microwave  region  of  the  electromagnetic  spec- 
trum (3  to  10  GHz).   A  distance  is  determined  by  observing 
the  time  needed  for  a  pulse  to  travel  from  a  master  tran- 
sponder located  aboard  the  survey  vessel  to  a  remote  tran- 
sponder en  shore  and  back  to  the  master  transponder. 
Knowing  the  average  velocity  of  the  electromagnetic  pulse, 
the  distance  D  is  then 

0  =  ^  (2.1) 

where  c  is  the  group  velocity  of  the  wave  packet  and  t  is 
the  two-way  travel  tine.   Short-range  systems  which  are  in 
wide  use  today  are  Bacal  Decca's  "Trisponder"  and  Motorola's 
"Mini-Banger."   These  systems  have  direct  range  readout  and 
are  readily  interfaced  into  a  navigational  computer  and  a 
data  acquisition  system.   Both  systems  are  active  and  user 
limited. 

Medium-range  positioning  systems  operate  in  the  1- 
to  5-MHz  frequency  range  of  the  electromagnetic  spectrum.   A 
distance  is  determined  by  measuring  the  phase  relationship 
between  transmitted  and  received  waves.   These  systems  are 
usually  referred  to  as  continuous  wave  systems  and  the 
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problem    cf   lane   ambiguity   must    be   addressed.      Ranges    are 
expressed  in  full  and  partial    lane   counts    where    a    lane    width 
w  is 

w  -i 

2  (2.2) 

where  A  is  the  wavelength  of  the  transmitting  frequency,  f, 
and  given  by 

A  «  f-  (2.3) 

Medium-range  systems  ccmmonly  in  use  today  are  Cubic 
Western's  "AEGO,"  Hasting  Raydist's  "Raydist,"  and  Odom 
Cff shore's  "Hydrotrack. " 

She  angle  of  intersection  associated  with  a  range- 
range  position  is  computed  frcm  a  simple  trigonometric  rela- 
tionship.  The  vessel's  position  P  (Fig.  2.1)  is  determined 
by  the  intersection  cf  the  ranges  from  the  left  and  right 
shore  stations,  R1  and  R2  respectively.   3  is  the  base  line 
distance  computed  between  the  two  known  shore  stations. 
Since  the  range  circles  from  the  shore  stations  intersect  at 
two  points,  it  is  necessary  for  the  plotter  to  recognize 
which  side  of  the  base  line  the  vessel  is  on  in  order  to 
eliminate  the  ambiguity.   The  angle  of  intersection  cf  the 
two  LOP's  (6)  is  given  by  the  law  of  cosines 

_  d2     pi  2     po2 

8  >  180°  -  Arc  cos  (  °  \ Kfo "jg       )  (2.4) 

In   qualitative    terms,   the   fix   is    strongest    when  8    approaches 
90°.      Host   hydrographic  specifications  limit   the    angle    cf 
intersection  from  a    linimum    of    30°    to  a    maximum   of    150°. 
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Figure  2. 1    Gecmetry  of  a  Range-Bange  Position 

2  •   Hyperbolic- Hyperbolic 

Hydrographic  positioning  by  hyperbolic-hyperbolic 
geometry  utilizes  the  intersection  of  two  hyperbolas  each 
generated  about  a  pair  of  shore  control  stations.   A  hyper- 
bola is  the  locus  of  points  in  which  the  difference  cf 
distance  from  two  fixed  points  is  always  constant.   A  three- 
station  hyperbolic  net  is  the  most  commonly  used  hypertclic 
mode  for  offshore  survey  (Fig.  2.2).   One  family  of  hyper- 
bolas (Red)  are  generated  about  a  master  station,  M,  and  a 
slave,  3;  while  a  second  family  of  hyperbolas  (Green)  are 
generated  with  respect  to  the  master  and  a  second  slave,  G. 
For  the  first  family  cf  hyperbolas,  the  control  points  K  and 
E  act  as  the  foci,  while  points  H  and  G  act  as  the  feci  fcr 
the  second  family. 
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Hyperbolic  location  methods  can  be  divided  into  two 
groups  based  on  the  electronic  principles  used  to  define  the 
distance  differences  [Ref.  7,  p.  87].   Loran  is  an  example 
of  a  pulse  system  in  which  the  differences  in  times  cf 
arrival  of  pulses  transmitted  ry  the  master-slave  combina- 
tions are  translated  into  distance  differences.   The  resul- 
tant position  has  no  lane  ambiguity  and  is  easily  resolved. 
The  second  method  of  hyperbolic  positioning  involves  meas- 
uring a  phase  difference  from  two  master-slave  combinations 
at  the  vessel's  position.   The  phase  difference  translates 
into  a  fractional  lane  count  which  in  itself  provides  an 
ambiguous  position.   This  ambiguity  is  resolved  by  using  a 
whole-lane  counter  which  is  initialized  at  a  known  geograph- 
ical point.   In  hyperbolic  positioning,  the  ship  is  in  a 
passive  mode  and  the  system  can  be  used  by  many  vessels. 

The  angle  of  intersection  between  the  two  hyperbolas 
can  he  computed  by  first  defining  the  following  guantities: 
S   is  the  length  of  red  base  line, 

S   is  the  length  of  green  tase  line, 
g 
R   is  the  distance  tetween  master  and  vessel's  position  ?, 

R  is  the  distance  from  red  slave  to  point  P, 

R  is  the  distance  from  green  slave  to  point  P, 

a  is  the  angle  between  lines  PM  and  PR,  and 

a  is  the  angle  between  lines  PM  and  PG. 

The  spacing  between  lanes  increases  with  distance 

from  the  master-slave  pair.   The  lane  widths  along  the  tase 

line  are 

w'  _  X         and        w'   X        ,_  _ 
r  "  ^L  9  =  9        (2.5) 

Then  the  lane  widths  at  any  point  P  are 

X        ,        Mj   ..  .  X, 


=  r  ,     1     \   and   w_  -    g  f   I  \ 

r  J-  {   sin  (a  72)  }         9   2  (  sin  (a/2)  }         (2.6) 
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where  tie  term  1/sin  (  a/2)  is  called  the  lane  expansion 
factor.   The  angle  oi  intersection  #g,  between  the  two 
hypertolas  is  then  given  by 


B  -  ar  *   °q 


(2.7) 


\  \       /  /   / 
\  \  \      /  /   / 

\  \  \    \  \         I    I   / 
\    \  /  / 


Figure  2.2   Geometry  of  a  Hyperholic-Hyperbolic  Position 


j. 


Ba nqe-A zimut h 


This  positiocing  geometry  is  used  for  nearshore, 
line-cf-sight  surveys.  One  LOE  is  generated  by  an  elec- 
tronic range  originating  from  a  transmitter  located  on  a 
shore  control  station.   A  microwave  system  is  commonly  used 
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in  this  arrangement  tut  systems  employing  a  laser  can  alsc 
be  used  for  short-range  work.   Another  LOP  is  generated  by 
fixing  an  azimuth  frcm  a  shore  control  station  to  the 
vessel.   A  second  control  station  is  used  for  an  initial 
azimuth  ty  the  observer.   Azimuth  determinations  can  be  made 
after  observing  directions  with  a  theodolite  as  an  observer 
tracks  the  moving  vessel. 

There  are  twc  ways  to  determine  a  range-azimuth 
position.   The  most  ccmmon  way  is  to  have  the  theodolite  and 
the  transmitter  occupy  the  same  shore  control  station. 
Hence,  the  angle  of  intersection,  g,  of  the  LOP's  is  always 
90°.   This  arrangement  is  commonly  used  by  the  National 
Ocean  Service  (NOS)  for  large-scale  nearshore  surveys. 

The  other  way  is  to  have  the  theodolite  and  the 
transmitter  occupy  two  different  control  points.   Then  the 
geometry  is  similar  to  that  of  the  range-range  position. 
The  angle  of  intersection,  3,  is  computed  by  trigonometric 
relationships  among  the  azimuth  of  a  line  between  the  shore 
stations,  the  observed  direction  to  the  vessel,  and  the 
measured  range  to  the  vessel. 

U .   Azimuth- Azimuth 

Azimuth-azimuth  positioning  geometry  is  used  for 
nearshore  high-accuracy  surveying.   Theodolites  are  set  over 
two  control  stations  en  shore.   The  vessel  is  sighted  on 
simultaneously  by  the  two  theodolite  observers,  generating 
two  visual  LOP's  whose  intersection  define  the  vessel's 
location.   Initial  azimuths  are  fixed  by  sighting  en  control 
stations  which  are  visible  to  the  observers. 

The  angle  of  intersection  for  an  azimut h-aziiruth 
position  is  dependent  on  the  geometric  relationships  tetween 
the  occupied  stations,  the  initial  stations,  and  vessel's 
position  (Fig.  2.3).   Assuming  that  theodolite  observers 
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occupy  stations  1  and  2,  and  initial  on  stations  3  and  'A, 
respectively,  the  observer  at  station  1  measures  angle  Y  and 
the  observer  at  station  2  measures  y    to  the  vessel.   The 

2 

angle  of  intersection,  8/  is  then  computed  by  first  deter- 
mining the  forward  azimuths,  measured  clockwise  from  the 
south,  ficm  stations  1  to  2  (a  ) ,  1  to  3  (a  ) ,  2  to  1 

12  13 

(  a  ) ,  and  2  to  4  (a  ) .   The  interior  angles,  3   and  9  ,  of 

21  2<»  12 

triangle    12P  are 

I  13  1  12 '  l*-°J 

and 

9     -   |a     +  y    -  a     I  ,~q, 

2  24  2  2  1  I* -yJ 

so   the    angle  of    intersection,    8    ,    at    the    vessel's    location 
is 

6  '  180°  "  <9,   +  V  (2.10, 


E.   C1ASEES  OF  EEBOBS 

All  hydrographic  positioning  measurements  are  subject  to 
error.   The  following  sections  discuss  categories  of  errors 
and  methods  used  to  treat  these  errors. 

1 .   Elunders 

Blunders  are  gross  mistakes  which  are  generally  due 
to  the  carelessness  of  the  observer.   Blunders  can  vary  in 
magnitude,  ranging  from  large  errors  which  are  easily 
detected,  to  small  errors  which  may  be  barely  distinguished. 
They  can  be  detected  by  making  repeated  observations  cr  by 
carefully  checking  the  data  in  the  processing  phase. 
Blunders  occur  in  various  forms  and  most  can  he  avoided  by 
carefully  planning  the  data  acguisition  process. 
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Figure  2.3    Geometry  of  an  Azimuth-Azimuth  Position 

Consider  the  following  as  an  example  of  a  blunder 
associated  with  range-range  geometry.   An  offshore  position 
is  to  he  determined  by  the  intersection  of  two  electronic 
LOP's  generated  from  transmitters  located  on  known  shore 
stations.   The  vessel  is  working  west  of  a  shoreline  that 
runs  generally  in  a  north-south  direction.   As  the  hydrc- 
grapher  faces  the  stations  from  sea,  the  southern  shore 
station  is  mistakenly  identified  as  left  and  the  northern 
shore  station  as  right.   The  resultant  offshore  position 
will  plot  to  the  east  of  the  base  line.   This  blunder  is 
readily  detected  and  can  be  easily  remedied. 

Not  all  types  of  blunders  are  so  easily  detected. 
Suppose  an  offshore  position  is  to  be  determined  ty  a 
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range-azimuth  fix.   A  range  and  an  azimuth  are  generated 
from  a  known  control  station  to  the  vessel's  position.   A 
second  control  station  is  used  to  fix  the  initial  azimuth;  a 
third  shore  control  station  is  located  10  meters  from  the 
initial  station  and  its  coordinates  are  mistakenly  used  for 
the  initial  station  in  plotting.   The  resultant  hydrographic 
position  is  in  error,  but  this  error  will  not  be  easily 
distinguished. 

Although  most  blunders  have  their  origin  in  human 
carelessness,  some  can  be  attributed  to  equipment  malfunc- 
tion.  For  example,  microwave  systems  which  generate  LOT'S 
are  known  to  become  unsteady  under  certain  conditions. 
Spurious  range  readings  resulting  from  signal  reflections 
can  te  recorded  as  true  positioning  data.   In  this  case,  the 
blunder  may  or  may  not  be  easily  detected. 

In  automated  data  acquisition  systems,  software  has 
teen  developed  to  detect  the  occurrence  of  anomalous  range 
readings.   By  inputting  a  course  and  speed  of  a  vessel  trav- 
eling along  a  line,  the  computer  can  determine  if  the 
recorded  position  is  valid  based  on  the  principle  of  dead 
reckoning.   If  the  recorded  position  is  found  to  be  invalid 
the  hydrographer  will  be  immediately  alerted  to  the  situ- 
ation and  can  take  action  to  remedy  the  problem.   In  non- 
automated  systems  the  principle  of  dead  reckoning  is  applied 
manually.   Given  the  course  and  speed  of  the  vessel,  the 
validity  of  the  position  can  be  checked  with  spacing 
dividers.   This  involves  checking  the  spacing  between  fixes 
recorded  before  and  after  the  position  in  question. 

Eefore  any  type  of  error  analysis  is  to  be  performed 
en  the  hydrographic  positioning  data,  it  is  essential  that 
all  blunders  be  identified  and  properly  treated.   In 
general,  careful  planning  coupled  with  thorough  checking 
will  micimize  the  occurrence  of  blunders. 
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2-   Systematic  Errors 

Systematic  errors  occur  with  the  same  sign,  usually 
of  similar  magnitude,  and  can  te  expressed  in  terms  of  a 
mathematical  model-   Systematic  errors  follow  a  defined 
pattern  and  occur  in  a  number  cf  consecutive  related  cbser- 
vations.   Eepetition  cf  measurements  does  nothing  tc  mini- 
mize their  effect.   Ir  the  case  of  hydrographic  positioning, 
systematic  errors  are  identified  and  modeled  by  calibration 
of  the  measuring  instrument  against  a  known  standard.   The 
following  is  a  brief  discussion  concerning  systematic  errors 
and  their  treatment  in  relation  to  hydrographic  positionirg 
eguipmer.t. 

a.   Theodolites 

In  nearshcre  surveys  the  theodolite  is  used 
primarily  for  range-azimuth  and  azimuth-azimuth  positicring. 
Systematic  errors  associated  with  the  theodolite  can  be 
classified  into  two  groups:  those  associated  with  the  phys- 
ical design  of  the  instrument  and  those  involving  the  geom- 
etry of  the  positioning  scheme.   Some  sources  of  systeiiatic 
errors  £Ref.  8]  associated  with  the  physical  characteristics 
cf  a  theodolite  are: 

i.   The  horizontal  circle  may  be  eccentric, 
ii.   Graduations  on  the  horizontal  circle  may  not  te 
uniform, 
iii.   The  horizontal  axis  of  the  telescope  (about  which 
it  rotates)  may  not  be  perpendicular  to  the 
vertical  axis  of  the  instrument, 
iv.   The  longitudinal  axis  of  the  telescope  may  net  te 
normal  to  the  horizontal  axis, 
v.   The  telescope  axis  and  the  axis  of  the  leveling 
bubble  may  net  be  parallel. 
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These  errors  are  usually  small  in  magnitude  and  can  be  elim- 
inated by    proper  adjustment  of  the  instrument  by  either  tie 
manufacturer  or  a  gualified  technician. 

The  field  hydrographer  has  ultimate  control  over 
the  geometric  systematic  errors  associated  with  a  theodo- 
lite.  In  range-azimuth  positioning  the  theodolite  and 
transmitter  may  occupy  the  same  horizontal  control  station. 
If  the  theodolite  is  not  set  directly  over  the  station  a 
resultant  systematic  error  will  occur  in  all  measurements. 
It  can  be  shown  that  these  errors  are  non- linear  but  do 
follow  a  mathematical  relationship.   Likewise,  if  the  trans- 
mitter is  not  located  directly  over  the  station,  a  similar 
type  cf  bias  occurs.   Depending  on  the  eccentricity  of  the 
theodolite,  the  vessel's  range  from  the  theodolite,  and  the 
scale  of  the  survey — these  errors  can  seriously  affect  the 
absolute  accuracy  of  the  offshore  positions. 

In  a  similar  fashion,  it  is  also  imperative  to 
position  the  target  directly  over  the  horizontal  control 
station  used  as  an  initial.   Failure  to  do  this  will  result 
in  an  error  which  will  be  propagated  to  offshore  positions. 

Many  situations  arise  in  the  field  where  it  is 
advantageous  to  set  a  transmitter  and  theodolire  over  a 
single  horizontal  control  station.   Freguently  it  is 
feasible  to  construct  a  platform  to  accommodate  both  instru- 
ments; in  a  case  where  it  is  not,  the  position  of  an  eccen- 
tric horizontal  control  station  near  the  original  station 
should  be  determined  and  that  station  used  for  the  location 
of  one  of  the  instruments.   The  theodolite  and  the  trans- 
mitter then  occupy  tvc  known  stations  and  the  geometric 
source  of  systematic  error  is  eliminated. 

b.   Electronic  Ranging  Systems 

The  systematic  errors  associated  with  electronic 
positioning  systems  are  complex  in  nature  and  functions  of 
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many  variables.   Munscn  [Hef.  9,  p.  4]  addresses  several 
problems  associated  with  short-range  systems  used  ir.  hydro- 
graphic  surveys.   The  most  common  problems  with  short-range 
systems  are  variation  in  range  and  calibration  drift  with 
time.   Variations  in  internal  eguipment  time  delays  in  the 
transmitter,  the  transponder,  or  the  receiver  can  induce 
errors  in  measured  ranges.   For  pulse  systems  such  varia- 
tions can  occur  due  to  temperature  dependence  of  components 
and  fluctuations  in  signal  strength  at  the  transponder. 
Multipath  effects  are  also  a  problem.   Under  some  circum- 
stances a  reflected  wave  and  the  directly  transmitted  wave 
arrive  with  a  phase  difference  of  180°.   Cancellation  or 
fading  of  the  directly  transmitted  signal  can  result. 

NOS  conducts  base  line  calibrations  of  short- 
range  positioning  systems  periodically  during  the  course  of 
a  survey  to  minimize  or  eliminate  systematic  error.   In  this 
process,  a  transmitter  and  receiver  are  each  placed  over 
control  stations  on  shore  and  the  measured  range  is  compared 
to  the  true  range.   In  this  way  the  systematic  error  is 
eliminated  by  zeroing  the  instrument  or  by  applying  a 
constant  correction  to  raw  data.   System  checks  are 
performed  daily  to  assure  there  is  no  drift  from  the  orig- 
inal calibration.   A  check  can  be  accomplished  by  comparing 
a  position  defined  by  the  ranging  system  to  a  known  fixed- 
point  position,  to  a  sextant  fix  position,  or  an  intersec- 
tion position. 

Munson  [ Bef .  9,  p.  5]  also  discusses  sources  of 
systematic  errors  associated  with  medium-range  systems.   Ihe 
most  significant  systematic  errors  occur  as  a  function  of 
position  due  to  varying  propagation  velocity.   The  medium- 
range  electronic  signal  propagation  velocity  depends  on  the 
surface  conductivity  and  transmission  path  (over  water,  over 
land,  or  over  different  types  of  land) .   Because  of  this 
dependence,  systematic  errors  as  a  function  of  position 
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cccur  at  different  effective  phase  velocities.   Knowing  the 
propagation  velocity  to  use,  or  the  phase  correction  tc  make 
as  a  furction  of  range,  is  a  problem.   Sky  wave  and  storm 
interference  also  pose  problems.   At  extreme  ranges  of  oper- 
ation, sky  wave  interference  can  affect  the  more  predictable 
ground  wave,  especially  during  nighttime  operations.   lane 
ambiguities  are  also  a  problem.   Most  systems  are  inherently 
ambiguous  and  must  be  zero  set  and  continually  monitored  for 
lane  jumps  or  loss  of  signal  which  results  in  the  loss  cf 
lane  count. 

NOS  uses  several  technigues  to  determine  the 
systematic  error  associated  with  medium-range  positioning 
systems.   These  technigues  involve  determining  a  whole  and 
partial  lane  count  for  phase  comparison  systems.   Two  of  the 
more  widely  used  techniques  are  comparison  of  three-point 
sextant  fix  positions  to  positions  determined  by  the  elec- 
tronic ranging  system  and  calibration  of  the  electronic 
system  at  a  fixed  point.   In  both  technigues  the  whole  lane 
counts  are  fixed  by  the  calibration;  correctors  to  the 
partial  lane  count  are  determined  and  applied  to  the  raw 
ranging  data. 

3 .   Ban dom  Errors 

Random  errors  are  chance  errors,  unpredictable  in 
magnitude  or  sign,  and  are  governed  by  the  laws  of  prob- 
ability [Ref-  10,  p.  1206  ]„   They  are  errors  which  remain 
after  blunders  and  systematic  errors  have  been  removed. 
Random  errors  result  from  accidental  and  unknown  combina- 
tions of  causes  and  are  beyond  the  control  of  the  observer. 
Greenwalt  [Ref.  12,  p.  2 ]  states  they  are  characterized  by: 
i.   Variation  in  sign;  positive  errors  occur  with 
equal  frequency  as  negative  ones, 
ii.   Small  errors  cccur  more  frequently  than  large  errors, 
iii.   Extremely  larce  errors  rarely  occur. 
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Kandom  errors  are  unique  to  specific  types  of  posi- 
tioning eguipment  and  vary  in  magnitude  depending  on  the 
precision  of  the  instruments  that  are  used.   The  following 
section  outlines  statistical  methods  for  their  treatment. 

C.   TEEA3BENT  OF  RANICH  EBBOBS. 

1 .   Cne-Di  mensional  Errors 

Certain  basic  statistical  quantities  must  first  be 
defined  in  the  analysis  of  random  errors.   Consider  a  vessel 
moored  securely  to  a  fixed  offshore  platform.   A  number  of 
ranges,  n,  from  a  microwave  transmitter  located  on  a  shore 
control  station  are  recorded.   The  mean  of  these  observa- 
tions is 

n  x. 

u  -  Z     — 
x   1.1  n  (2.  11) 

where  x   represents  an  individual  observation.   The  standard 
error,  s,  of  the  observations  is  then 


S"M^A  (xi-u*)2  (2-12) 


where  the  quantity  (x  -  u  )  is  referred  to  as  the  residual, 

i       x 

cr   true    error,    v    ,    of  a   particular   observation.      As   n    gets 

very    large,    the    factor    1/n  can    be   substituted  for    1/(n-1)    in 

Equation   2.12.       likewise,    in    treating    the    large  sample,  a 

can    be    substituted   for  s   and   h   for  V    $    where  y    and    g     are 

x 

the  mean  and  standard  error  of  the  entire  population. 

It  is  of  interest  to  determine  the  probability  of 
occurrence  of  a  particular  observation.   The  normal  cr 
Gaussian  distribution  equation  relates  the  residual  of  a 
particular  random  variable  with  the  probability  of  its 
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occurrence,    and    is   given   by 


P(v)  . 


-<£* 


(2.13) 


The    plot   of    this   equation   yields    the   normal   distrifcuticn 
curve    (Fig.    2.4).       The   height    cf    the   curve    above    the 
vertical  axis    is   proportional    to   the   probability    of    a 
particular   error   occurring. 

The   probability   of   a  residual   falling   between    asy 
two    residuals    v      and    v      can    be   computed    by    integrating 

1  2 

Equation   2.13    as 


P(v)  = 


-<& 


>ffir 


dv 


(2.14) 


Figure    2.4        The    Normal    Distribution 


This   integral   is   difficult    to   evaluate    analytically 
so    tables    have    been    compiled    to    aid    in    computations.       Fcr 
v      =    *a    and    v      =   -a,    it   can   be    shown    that    P(v)    =    0.6827.      In 

1  2 

ether    words,    the    probability    that    a    particular    observation 
will    fall   within    +    1a      of    the    mean   is    68.27    percent. 
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Eeturning  to  the  example  of  the  vessel  moored  tc  the 
offshore  platform,  th€  mean  and  the  standard  errcr  fcr  the 
observations  are  easily  computed.   With  this  information  and 
Equation  2.14,  the  probability  of  a  range  error  falling 
within  specified  limits  can  be  computed.   Conversely,  by 
fixing  a  probability,  the  associated  limits  of  the  range 
error  can  he  computed  .   In  statistical  terms,  a  particular 
observation  will  fall  within  specified  limits  with  a  certain 
confidence. 

Actual  values  of  one-dimensional  standard  errors  for 
hydrographic  positioiing  equipment  are  a  subject  of  debate 
betweer  manufacturers  and  users.   Some  manufacturers  of 
microwave  positioning  equipment  claim  standard  errors  of  +1 
meter.   Cn  the  other  hand,  Munson  [ Ref .  9,  p.  6]  states  that 
microwave  systems  demonstrate  accuracies  of  3  meters  at 
short  ranges  but  show  larger  errors  at  ranges  of  15  km  and 
greater.   NOS  assumes  a  3-meter  standard  error  in  all  of  its 
short-range  accuracy  computations.   It  is  apparent  that 
further  study  is  needed  to  adequately  define  the  nature  of 
errors  associated  with  electronic  positioning  equipment. 

Waltz  [Ref.  13]  performed  an  extensive  study  to 
detemine  the  pointing  error  of  a  Wild  T-2  theodolite.   His 
results  showed  that  the  pointing  error  associated  with  this 
instrument  under  hydrographic  survey  conditions  was  about 
1.3  meters  and  was  independent  of  distance. 

2 .   Iwo-Di mensional  Errors 

The  intent  of  this  paper  is  to  apply  statistical 
methods  developed  by  ethers  tc  a  hydrographic  data  set 
containing  two-dimensional  errors  which  are  defined  by  two 
random  variables.   lengthly  and  complex  derivations  are  not 
presented.   Burt  [Eef.  3]  and  Heinzen  [Ref.  2]  show  adeguate 
derivations  of  formulas  associated  with  two-dimensional 
errors  ar.d  can  be  referenced  for  full  details. 
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Ihe  following  assumptions  are  made  concerning  two- 
dimensional  errors  associated  with  intersecting  LOP's: 

i.   The  random  encrs  of  each  LOP  are  normally 
distributed, 
ii.   Systematic  or  bias  errors  have  been  removed  from 
the  observations, 
iii.   The  intersecting  LOP's  are  coplanar. 
iv.   The  error  LOP's  are  parallel  to  the  exact  LOE's. 

In  developing  a  usable  mathematical  model  for  accuracy 
determinations,  the  four  assumptions  hold  to  a  high  degree 
for  all  hydrographic  positioning  geometries. 

Consider  again  the  vessel  moored  to  a  fixed  offshore 
platform.   Assume  two  ranges  are  measured  from  two  different 
shore  control  stations  at  the  same  time  and  that  the  range 
readings  are  uncorrelated.   The  observation  of  this  pair  of 
ranges  is  repeated  many  times.   After  a  large  number  of 
observations,  the  means  and  standard  errors  of  the  indi- 
vidual ranges  are  determined.   Suppose  the  mean  ranges,  or 
the  actual  LOP's,  intersect  at  an  angle  of  90°  and  that  the 

computed  standard  errors  are  egual  (a  -a     )•   If  each  data 

1    2 
pair  (x  ,y  )  is  plotted,  the  spread  of  points  about  the  mean 

coordinates  results  in  a  circular  cluster  (Fig.  2.5).   A 

higher  density  of  points  occurs  near  the  intersection  of  the 

mean  ranges  and  the  density  of  points  decreases  outward  from 

the  intersection  of  the  mean  ranges. 

In  this  special  case,  which  is  called  a  circular 

normal  distribution,  the  probability  of  a  point  falling 

within  a  specified  radius,  R,  from  the  intersection  of  the 

mean  ranges  is 

l^TJ  (2.15) 

P(R)  »  1  -  e   c 
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where  a  =  a  =  °~     and  is  defined  as  the  circular  standard 

1       2      C 

error.   Osing  Eguaticn  2.15,  R  can  be  computed  by  fixing 

P (R) ,  cr  ccnversly,  P (R)  can  be  computed  by  fixing  R. 

letting  R  =  a   =  a   =  a  ,    then  P(R)  =  0.3935.   In  other 

1     2     c 
words,  39.35  percent  cf  all  errors  in  a  circular  normal 

distribution  are  not  expected  to  exceed  the  circular  stan- 
dard error  [Ref.  12,  pp.  25t26]. 


TO  LEFT  STATION 


LOP  2 

TO  RIGHT  STATION 


Figure  2.5   Circular  Normal  Distribution 


In  the  case  where  the  two  uncorrelated  LCPfs  inter- 
sect at  an  angle  other  than  90°  or  a  ^ a  ,  the  contours  cf 

1      2 
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equal  density  are  ellipses  centered  about  the  point  defined 
by  the  intersecting  ICP's  (Fig.  2-6).   The  two-dimensional 
probability  density  function  becomes  [Ref.  1,  p.  136] 


P(Vvy)=^e 


K2 

7 


(2.16) 


Figure  2-6   Error  Ellipse  Formed  by  Two  Oncorrelated  LCP's 


where 

v   is  the  residual  in  the  direction  of  the  semi-major  axis 
x 


of  the  error  ellipse, 


v   is  the  residual  in  the  direction  of  the  semi-minor 
axis, 
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a     is  the  standard  error  in  the  direction  of  the  semi- 
major  axis, 

a   is  the  standard  error  in  the  direction  ox  the  semi- 
y 
miner  axis, 

and 

vx   y2 

k2  =  4+4  (2*1?) 

The  solution  of  Equation  2.16  with  values  of  K  for  different 

P*s  yields  the  results  in  Table  I  [Ref.  12,  p.  23].   For  a 

39.35  percent  probability,  the  axes  of  the  ellipse  are 

1.0000  a   and  1.0000  a  ;  for  a  50  percent  probability,  the 

axes  are  1.1774 a   and  1.1774  a  - 

x  y 


TABLE  I 
Values  of  the  Constant  X 

PROBABILITY  K 

39.35%       j  1.0000 

50.00%  1.1774 

63.21%  1.4142 

90.00%       f  2.1460 

99.00%       |  3.0349 

99.78%       I  3.5000 


The  error  ellipse  can  be  used  for  accuracy  computa- 
tions by  developing  relationships  for  a   and  a   in  terms  of 

x  y 

the    initial  information  a    ,    a   ,    and     6  .      Bowditch    [Bef.    10, 

1  2 

p.    1213]  gives    the   following   equations  for   independent    iOP's 
relating   these    quantities: 

°x  =  2iik    {^♦a»Mlo{+aJ)i-4*^6oV   }       (2.18) 
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and 


^Sin  D       1      2  12  12        U«  iyj 


'y- 


In  these  equations,  8  is  assumed  to  be  the  acute  angle 
between  the  LOP's. 

In  certain  special  cases,  the  above  equations  take 
on  more  manageable  forms.   In  range-range  and  azimuth- 
azimuth  positioning  it  is  often  assumed  that  a      ~  a      =  a. 

1       2 

Equations  2.18  and  2.19  then  reduce  to 

a     -      ■  /T   B 

°x  "  2sin(>5B)  °  (2-20) 


and 


ay  "     2cosTW  (2.21) 

In    the   concentric   range-azimuth   case,    a     *    a     $    and  8 

1  2 

equals    90°.      Equations   2.18    and    2.19    then   simplify    to 

a       a 
x  -^  (2.22) 


and 


a     -  o 
y         2  (2.23) 


where  a      >  a      and    a     >  a    • 

12  X  y 

The  case   for   correlated   LOP • s  is    more  complex.      The 

calculation   of    0"     and  a      involves   a    coordinate  transf oria- 

x       y 

tion  from  a  linear  skewed  coordinate  system  to  an  uncorre- 
cted rectangular  cocrdinate  system.   The  following 
discussion  is  taken  from  Heinzen  [Ref.  2,  pp.  49-53]. 

Assume  a  hydrcgraphic  position  is  established  by  the 
intersection  of  two  correlated  LOP's  (Fig.  2.7a).   LOP  1  and 
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LOP  1 


LOP  2 


LINEAR 

SKEWED 

A 


LOP  2        X 


CORRELATED 

RECTANGULAR 

B 


UNCORRELATED 
RECTANGULAR 
C 


Figure  2.7   Coordinate  Transformations  for  Correlated  LOP's 

LOP  2  are  the  coordinate  axes  in  the  skewed  coordinate 
systen,  with  standard  errors  a   and  a  .   The  semi-major  arid 

1         2 

semi-minor  axes  of  tl€  error  ellipse  are  not  coincident  with 
the  skewed  coordinate  system  axes.   The  correlation  coeffi- 
cient between  the  two  LOP's  is  p  .   Assume  g-   >  a  • 

12  i       2 

Ihe  standard  errors  and  correlation  coefficient  in  a 

correlated  rectangular  coordinate  system  with  axes  A  and  B 
must  now  be  determined.   A  coordinate  transformation  from 
the  skewed  system  to  the  correlated  rectangular  system  nust 
be  made  yielding  the  standard  errors  along  the  new  coordi- 
nate axes  (Fig.  2.7b) 


°l   s  c-inZg  (a*  +  2p  a  a  cos8  +  a2)  -  a2 
a   sin  b    i     12  i  2       2     2 


(2.24) 
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and 


ab  -  a2 


(2.25) 


The  correlation  coefficient  in  the  correlated  rectangular 
system  is 

o  a  a     2  k 

Pflh  s  (  r2-  cosS  +  P     )   {   1  +  p    T-^cosB  +  (-^-)cos2B>"^     ,_    ,,. 
aD  a  12  12  o  0  (2.  2d) 

To    detemine    a      and  <J_  ,    a   second   coordinate    transformation 
x  y 

lust    be   performed    from   the   correlated   rectangular   system    to 
an   uncorrelated   rectangular   system   with   axes    X   and    Y    (Fig. 
2.7c).      The   semi-major   and   semi-minor  axes    of   the   error 
ellipse   are   then 


3T^T   /,     /.    ggTEZHb? 


l  +  ,/l r-TJ m-r  (2.27) 


i  °i  ♦  °b  ' 


and 


v        Y    a        b        x 


(2.28) 


When  P,   =  0.  these  equations  become  identical  to  the 

12  3 

simplified    versions    in   Bowditch   £Ref.    10]. 

The    orientation   of    the    semi-major   and   semi-miner 
axes    relative    to    the    intersecting   LOP's    is    the   third    param- 
eter   which    fixes    the    error   ellipse.       The   angle  9     (Figs.    2.6 
and    2.7)    is   measured   counter-clockwise   from   LOP    1     to    the 
semi-najcr   axis    of    the  error   ellipse  [Hef.     11]   and    is    given 

by 

.  o2sin(2B)  +  2o     a  a  s1n(S) 

9  *  i  arctan  {  -1 ■ ^-^ > 

L  a2cos(2B)  +  2p    a  a  cos(6)  +  a2                   (2.^9) 

1  12     1     2  2 

For    the   special   case    cf    a     =    a     and     p      =0, 

1  2  12 

9  *  I  (2.20) 
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The  orientation  of  the  error  ellipse  in  an  orthogonal  coor- 
dinate system  can  be  represented  by  adding  or  subtracting  9 
to  the  orientation  of  LOP  1.  Care  must  be  taken  on  deter- 
mining the  quadrant  of  the  outcome.  As  a  general  rule,  the 
error  ellipse  always  lies  within  the  acute  angles  termed  hy 
the  intersecting  LOP's. 

The  orientation  and  dimensions  of  the  error  ellipse 
provide  a  useful  index  for  evaluating  the  accuracy  of  a 
hydrographic  position.   Its  greatest  attribute  is  that  it 
accurately  represents  the  error  distribution  about  the 
intersection  of  two  ICP's  in  terms  of  a  fixed  probability. 
It  is  interesting  to  examine  the  variation  in  the  relative 
dimensions  and  orientations  of  error  ellipses  as  they  vary 

in  a  range-range  configuration  with  a     =  a   =  a   (fig-  2.8) . 

1     2 
The  dimensions  of  the  ellipses  are  specified  by  Equations 

2.20  and  2.21  and  a   and  o   are  functions  of  g  only  for 

x       y 
fixed  a .   Therefore,  the  dimensions  of  the  ellipses  remain 

constant  along  a  contour  of  constant  g  ;  only  the 

orientation  changes.   A  line  of  constant  g  is  a  circle 

which  includes  stations  L  and  5.   Note  that  the  dimensions 

of  the  ellipses  for  8  * s  of  30°  and  150°  are  identical.   Ihe 

ellipses  about  the  90°  angle  of  intersection  contour  are 

circles  and  represent  the  strongest  possible  positions  in 

this  scheme.   With  varying  8  fs,  the  directional  nature  of 

the  distribution  can  be  noted. 

3  •   Circular  Precision  Indexes 

Although  the  error  ellipse  gives  a  true  representa- 
tion of  the  error  distribution  about  a  hydrographic  posi- 
tion, its  use  has  certain  drawbacks.   The  characteristics  of 

the  ellipse  must  be  specified  by  the  three  quantities  a  , 

x 
a  ,  and  9.   A  single  figure  for  evaluating  the  positional 

accuracy  cannot  be  used.   Greenwalt  [Ref.  12,  p.  26]  states 

that  when  a      and  0  are  not  equal,  a  circular  error 

X  y 
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Figure   2-8        Error    Ellipses    Around   a    Range-Range    System 
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distribution  can  be  substituted  for  the  elliptical  distribu- 
tion.  This  substitution  can  be  satisfactory  for  error 
analysis  within  certain  ay^ax  ratios.   However,  when  this 
ratio  is  small  the  distortion  introduced  by  the  circular 
distribution  may  beccie  misleading. 

a.   Root  Mean  Square  Error 

The  terms  radial  error,  root  mean  square  error, 

and  d     are  identical  in  meaning  when  applied  to  t wo- 
rms 

dimensional  errors  [ Bef .  10,  p.  1229]-   The  term  d     is 

L  r  J  rms 

defined  as  the  sguare  root  of  the  sum  of  the  squares  of  the 
standard  errors  along  the  major  and  minor  axes  of  the  error 
ellipse.   That  is 

d   -  J  a*  +  a2 
rms  *  x    y 

(2.21) 

where   a      and    a     are    given   by    Equations   2.18    and   2.19. 
x  y 

A  more   direct   form    of  2.31    is    given    by   [Ref.    2,    p.    54] 


d         =       — --^7  l/o2   +  a2  (2.22) 

for    uncorrelated   LOP's,      For   range-range  and   azimuth-azimuth 

positioning,    with    cr     =   a      =    a    ,    Equation  2.32    reduces    to 

1  2 

drms  *     slTifi  °  (2-33) 

For    range-azimuth   positioning,     8  =    90°   and    Equation    2.32 
becomes 


™s  1         2 


(2. 34) 


The   mere   general  fori   of   Equation   2.32   for    both    correlated 
and    uncorrelated   LOP's   [fief.    2,    p.    59]  is 


d™.s  ■  %m  \ft*  °! +  2°12v2cos8         (2--5> 


where     p       is   the   correlation   coefficient. 
12 
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An    error    circle    with    a   radius    of    one    d  can    be 

rms 

constructed  about  the  intersecting  LOP's  (Fig.  2.9).   Two 

d     is  the  radius  cf  the  error  circle  obtained  using  two 
rms 

times  the  values  of  a     and  a   in  Equation  2.31.   For  an 

x  y 

elliptical   error  distribution,    the    probability  associated 

with    a    specific    value   of   d  varies   as   a    function    of    the 

rms 


eccentricity  of  the  error  ellipse  (Table  II) .   The  prob- 
ability associated  with  one  d     varies  from  63.2  percent 

rms 

68.3  percent,  while  the  probability  associated  with  twc  d 
varies  between  95.4  percent  and  98.2  percent. 


Figure  2.9   The  d    Error  Circle 
s  rms 


NOS  uses  d_„,  as  an  accuracy  specification. 

rms  ■*       fc 

Umbach   [Eef.    14,    p.    4-25]    states    that    super    high    freguency 
direct    distance    measuring   systems    would    be    used    only    when 
the    value   cf  dTms     is    less   than   or   egual    to: 
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i.   0.5  mm  at  the  scale  of  the  survey  for  scales  of 

1:20,000  and  smaller, 

ii.   1.0  mm  at  the  scale  of  the  survey  for  1:10,000 

scale  surveys,  or 

iii-   1.5  mm  at  the  scale  of  the  survey  for  scales  of 

1:5,000  and  larger. 

The  major  advantage  of  using  d    as  a  precision 

rms 

index  is  its  ease  of  computation.   Some  hydrographers  draw 

analogy  between  the  varying  probability  associated  with  ore 

d     (63.2  percent  to  68.3  percent)  and  the  fixed  prct- 
rms 

ability  associated  with  a  one-dimensional  standard  error 


In  fact,  d     has  very  little  statistical 
rms 


(68.3  percent) 

meaning.   The  obvious  problem  with  using  dras  as  a  precision 
index  is  the  varying  trobability  associated  with  the  error 
circle-   For  this  reason  Greenwalt  [Ref.  12,  p.  31]  recom- 
mends against  its  use. 


TABLE    II 

Probabilities  Associated   Hith    d 

rms 

PROBABILITY 

0              a            LENGTH    OF 
y              x                       rms                       rms                     rms 

0.0 

1.0 

1.000         | 

0.683 

0.954 

0.1 

1.0 

1.005 

0.682 

0.955 

0.2 

1.0 

1.020 

0.682 

0.9  57 

0.3 

1.0 

1.042 

0.676 

0.961 

0.4 

1.0 

1.077 

0.671 

0.966 

0.5 

1.0 

1.  118 

0.662 

0.969 

0.6 

1.0 

1.  166 

0.650 

0.973 

0.7 

1.0 

1.220 

0.641 

0.977 

0.8 

1.0 

1.280         j 

0.635 

0.980 

0.9 

1.0 

1.345 

0.632 

0.  981 

1.0    j 

1.0 

1.414 

0.632 

0.982 
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b.   Circles  cf  Equivalent  Probability 

Burt  [Bef.  3]  presents  a  method  for  translating 
ellipses  of  equivalent  probability  into  circles  of  equiva- 
lent probability.   Tc  utilize  this  method,  it  is  first 
necessary  to  compute  the  eccentricity  of  the  error  ellipse, 
c,  by  the  equation 

a 

c  --* 

ax  (2.36) 

where  a    >  a    . 
x    y 

Harter  [Bef.  15]  compiled  Tables  III  and  IV 
which  are  taken  from  Bowditch  [Ref.  10,  p.  1215]-   Harter's 
data  are  given  in  terns  of  the  eccentricity,  c,  a  parameter, 
K,  and  a  probability,  P.   The  parameter,  K,  when  multiplied 
by  a     gives  the  value  of  the  radius,  R,  of  the  circle  of  the 
corresponding  probability  shown  in  Table  III.   That  is, 

R  s  K  °X  (2.37) 

The  probability  of  a  point  falling  inside  a 
circle  cf  specified  radius  can  be  computed  by  entering  Table 
III  with  c  and  K  as  arguments.   Given  a  fixed  probability,  K 
is  deternined  by  entering  Table  IV  using  c  and  P  as  argu- 
ments.  Ihe  radius  of  the  probability  circle  is  then 
computed  using  Equation  2.37. 

Using  confidence  ellipses  has  certain  advantages 
ever  confidence  circles  of  equal  probability.   First,  the 
directional  nature  of  the  true  error  distribution  is  net 
represented  in  the  confidence  circle  method  even  though  both 
methods  give  an  accurate  measure  of  confidence.   Second,  the 
area  of  the  confidence  ellipse  is  always  less  than  cr  equal 
to  the  area  of  the  confidence  circle.   The  area  of  a 
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TABLE  III 
Probabilities,  Given  c  and  H 


K    \ 

a  0 

0.  1 

0.  2 

a  3 

0.  4 

0.5 

1 

0.  6 

0.  7 

0.8 

a  9 

1.  0 

0.  1 

.  079C557 

.  0443987 

.  0242119 

.  0164176 

.  0123875 

.  0099377 

.  0082940 

1 
.  0071157 

.  0062299 

.  0055400 

!   0049875 

0.2 

.  1585194 

. 1339783 

.  0884533 

.  0628390 

.  0482413 

.  0390193 

.  0327123 

1  .  0281415 

.  0246824 

.  0219757 

1  .  0198013 

0.  3 

.  2358228 

.  2213804 

.  1739300 

.  1318281 

.  1039193 

.  0651535 

.  0719102 

1  .  0621336 

0546598 

.  0487639 

0440025 

0.  4 

3108435 

.  3010228 

.  2635181 

.  2139084 

.  1742045 

.  1451308 

.  1237982 

!  .  1076237 

j  .  0950495 

.  0850320 

.  07C8837 

0.  5 

3829249 

.  3755884 

.  3481790 

.  3003001 

.  2532953 

.  2152880 

.  1357448 

|  .  1020829 

.  1443941 

.  1290236 

.  1175031 

0.  6 

4514938 

.  4457708 

.  4255605 

.  3S46374 

.  3357384 

.  2914682 

.  2548177 

.2251114 

2009797 

.  1811783 

1 047293 

0.7 

.  5160727 

.  5115048 

|  .  4900683 

.  4633258 

.  4170862 

.  3699305 

.  3280302 

.  2925654 

2P29373 

I  .  2381583 

.  2172955 

0.8 

.  5762892 

.  5725957 

.  5604457 

.  5349387 

.  4941882 

.  4474207 

|  .  402562S 

.  3627122 

.  3283453 

1  .  2939700 

!  .2733510 

0.  9 

6318797 

.  6288721 

.  6191354 

.  5993140 

.  5651564 

.  5213998 

.  4759375 

. 4333628 

.  3953279 

.  3C20135 

1  .  3330232 

1.0 

.  6826895 

.  6802325 

.  6723580 

. 6568242 

.  6291249 

.  5900953 

.  5461319 

1  .  5025790 

.  4621421 

4257553 

j  .  3934693 

1.  1 

.  7286679 

.  7266597 

.  7202682 

.  7079681 

.  6859367 

.  6324489 

.  6116316 

.  5687467 

.  5272462 

.  48S7S73 

.  4539256 

1.2 

.  7698607 

.  7682215 

.  7630305 

.  7532175 

.  7359558 

.  7079973 

.  6714269 

.  6306168 

.  5893494 

.  5498736 

I  .  5132477 

1.3 

.  8063990 

.  8050648 

.  8008554 

.  7929968 

.  7793550 

.  7567265 

.  7249673 

1  .  6873122 

6474394 

.  6079822 

[  .  5704426 

1.  4 

.  83S4867 

.  8374049 

.  8340013 

.  8277048 

.  8169851 

.  7989288 

.  7720S89 

.  7383089 

.  7007900 

.  6623035 

j   6216339 

1.  5 

.  8663856 

.  8655127 

.  8627728 

.  8577362 

.  8493071 

.  8350816 

.  81292S7 

.  7833962 

.  7489500 

.  7122546 

I  .  6753475 

1.  8 

.  8904014 

.  8897008 

.  8875060 

.  8834914 

.  8768644 

.  8657559 

.  8478393 

'  .  3226246 

.  7917194 

.  7574708 

.  7219627 

1.  7 

.  9108691 

.  9103102 

.  9085619 

.  9053766 

.  9001740 

.  8915536 

.  8773116 

.  3562471 

.  3291137 

1  .  7977882 

74>-2S39 

1.  3 

.  9281394 

9276964 

9263125 

.  9237989 

.  9197275 

9130680 

.  9019110 

'.  3846624 

.  8613238 

1  .  3332175 

.  3021013 

1.  9 

94256C9 

.  9422182 

.  9411299 

9391586 

.  9359855 

.  9308615 

.  9222277 

9033609 

.  3886731 

i  8639149 

.  8355255 

2.  0 

.  9544997 

.  9542272 

9533775 

9518415 

.  9493815 

.  9454546 

93S8418 

.  9278799 

.  9U5762 

.  5901495 

.  8046647 

2.  1 

.  9642712 

.  9640598 

.  9634011 

.  9622127 

.  9603170 

.  9573205 

.  9522999 

.  9437668 

9305013 

1  .  9122714 

.  8H97495 

2.2 

.  9721931 

.  9720304 

.  9715237 

9706109 

.  9691597 

.  9668845 

.  9631017 

.  9565522 

.  0459386 

,  .  9306821 

.  9110784 

2.  3 

.  9785518 

.  9784275 

.  9780408 

9773450 

.  9762419 

9745239 

.  9716934 

9667306 

.  9533739 

9458P35 

.  9289946 

2.  4 

.  9336049 

9835108 

9832130 

.  9826918 

.  9818594 

9805703 

.  9784661 

.  0747495 

.  9632698 

.  9580804 

.  942S652 

2.  5 

.  9875807 

.  9875100 

9872900 

.  9868953 

.  9862720 

9853112 

9337569 

0810035 

.  9760522 

.  9670 136 

.  9560631 

2.8 

.  9906776 

.  9906249 

.  9904612 

.  9901674 

.  9897045 

.  9889934 

.  9678527 

.  9858331 

.  9821023 

.  9756069 

.  9659525 

2.  7 

9930661 

9930271 

9929062 

.  9926894 

.  9923483 

.  9918260 

.  9909944 

.  9S9526S 

.  9867530 

.  9817837 

.  97387S6 

2.8 

.  9948S97 

.9948612 

.  9947727 

-  9946141 

.  9943649 

.  9939842 

.  9933821 

9923249 

.  9902888 

.  9864676 

. 0801589 

2.9 

.  9962684 

.  9962477 

.  9961834 

.  9960684 

.  9958878 

.  9956126 

.  9951798 

.  9944246 

.  9929432 

9900803 

.  9«50792 

3.  0 

.  9973002 

.  9972853 

.  9972391 

.  9971564 

.  9970266 

.  9968294 

.  9965205 

.  9959854 

.  9949274 

.  9927925 

.  9333910 

3.  1 

.  9980648 

.  9980542 

.  9980212 

.  9979622 

.  9978699 

.  9977296 

9975109 

.  9971348 

9963851 

.  9048168 

.  9913113 

3.2 

.  9986257 

.  99861S2 

.  9985949 

.  9985533 

.  9984880 

.  9983892 

.  9982356 

.  9079733 

.  9974478 

.  9063105 

.  9040240 

3.  3 

.  9990332 

.  9990279 

.  9990116 

.  9989824 

.  9989368 

.  9988677 

.  9987607 

.  9965792 

.  9932147 

.  9974004 

.  9956322 

3.  4 

.  9993261 

.  9993225 

.  9993112 

.  9992909 

9992593 

9992115 

.  9991376 

.  9990129 

.  9987626 

.  9981S68 

.  9060113 

3.5 

.  9995347 

.  9995323 

9995245 

.  9995105 

9994888 

9994559 

.  9994053 

9993204 

.  9991502 

.  99S7480 

.  0073125 

3.6 

.  9996818 

.  9996801 

.  9996748 

.  9996653 

9996505 

.  9996231 

.  9905933 

.  9995364 

.  9994218 

.  9901442 

.  0984662 

3.  7 

.  9997844 

,  9997832 

.  9997797 

.  9997733 

.  9997633 

.  9997482 

.  9997251 

.  9996367 

.  9990102 

.  9994208 

.  9989352 

3.  8 

.  9998553 

.  9908545 

.  9993522 

.  9998478 

9998412 

.  9998311 

.  9998157 

-  0997902 

.  9997396 

.  9996119 

.  9902682 

3.  9 

.  9999038 

.  9999033 

.  9999018 

.  9998989 

9998945 

.  9998378 

.  9998776 

.  9998606 

.  9908276 

.  9997426 

.  0095020 

4.  0 

9999367 

9999363 

.  9999353 

.  9999334 

9999305 

.  9999261 

.  9999195 

. 9999085 

.  9998870 

.  9998309 

.  0006645 

4.  1 

.  9999587 

.  9999585 

.  9999578 

.  9999566 

9999547 

.  9999519 

9999475 

.  9999404 

.  9999266 

.  9998900 

.  9997763 

4.2 

.  9999733 

.  9999732 

.  9999727 

.  9999720 

.  9999707 

.  9999689 

.  9009661 

9999616 

.  9999527 

.  9099292 

.  9993523 

4.  3 

.  9999829 

.  9999828 

.  9999826 

.  9999821 

.  9999813 

.  9999801 

.  9099733 

.  9099754 

.  9999698 

.  9999548 

.  9090034 

4.  4 

.  9999892 

.  9999891 

.  9999889 

.  9999S86 

.  9999881 

.  9999874 

0999863 

.  9999845 

.  9999809 

.  9900715 

.  9999375 

4.  5 

.  9999932 

.  9999932 

.  9999931 

.  9999929 

.  9999925 

.  9999921 

.  9999914 

.  9y99902 

.  9999881 

.  9999822 

.  9999599 

4.  6 

.  9999958 

.  9999957 

.  9999957 

.  9999955 

.  9999954 

.  9999951 

.  9000947 

.  9999939 

.  9990926 

.  9909389 

.  9999746 

4.  7 

.  9999974 

.  9999974 

.  9909973 

.  9999973 

.  9999971 

.  9999970 

.  9990067 

.  9009963 

.  9999955 

.  9909932 

.  9999840 

4.  8 

.  9999984 

.  9999984 

. 9999984 

.  9999983 

.  9999983 

.  9999982 

.  9999930 

.  0999977 

9909972 

.  9999950 

.  9999901 

4.  9 

.  9999990 

.  9999990 

. 9999990 

.  9999990 

.  9999990 

.  9999989 

.  90999S8 

.  9909986 

.  9999983 

.  9999075 

.  0000039 

5.  0 

.  9999994 

.  9999994 

9999994 

.  9999994 

.  9999994 

.  9990993 

.  0999993 

.  0999992 

.  9099990 

.  0909935 

.  0000063 

5.  1 

.  9999997 

.  9999997 

.  9999997 

.  9999996 

9999996 

.  9999996 

.  9999996 

.  9000995 

.  9999904 

.  0990991 

,  9999978 

5.  2 

.  9999998 

.  9999998 

.  9999998 

.  9999998 

.  9999998 

.  9999998 

.  9999998 

.  9009997 

.  0999997 

.  9999995 

.  9099987 

5.  3 

.  9999999 

.  9999999 

.  0999999 

.  9999999 

.  9999999 

.  9999999 

.  9999999 

.  9999998 

.  0999998 

.  9900097 

.  0090092 

5.  4 

.  9999999 

.  9999999 

.  9999999 

.  9999999 

9999999 

.  0999999 

.  9999999 

.  9909999 

.  9099999 

.  9909998 

.  9999995 

5.  5 

1.  OOOOOOO 

1.  0000000 

1.  0000000 

1.  0000000 

1.  0000000 

1.  0000000 

1  0000000 

.  9999999 

.  9909999 

.  9909999 

.  9909997 

5.  6 

1.  0000000 

I.  0000000 

.  9999999 

.  0090098 

5.  7 

I.  OOOOOOO 

.  9900999 

5.  8 
5.9 

6.  0 

i 

1  OOOOOOO 
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TABLE  17 
Radii  cf  Circles  Given  c  and  P 


^V     c 

p    >. 

0.  0 

0.  1 

0,2 

a  3 

0.  4 

0.  5 

0.  6 

0.  7 

0.  8 

0.  9 

!  L0 

.5000 
.7500 
.9000 

a  67449 
1.  15035 
1.  644S5 

0.  68199 

1.  15473 
1.  64791 

0.  70585 

1.  16825 
1.  65731 

0.  74993 

1.  19246 
1.  67383 

0.  80785 

1.  23100 
1.  69918 

0. 

1: 

87042 
2S534 
73708 

0.  933C5 

1.  35143 
1.  79152 

0.  99621 

1.  42471 
1.  86253 

1.  05769 
1.  50231 
1.  94761 

I  1.  11807 

,  1.  58271 

2.  04236 

1.  17741 
I  1.66511 

2.  14597 

.9500 
.9750 
.9900 

1.  95996 
Z  24140 
Z  57583 

1.  96253 

2.  24365 
2.  57778 

1.  97041 

2.  25053 
2.  58377 

1.  98420 

2.  26255 
2.  59421 

2.  00514 
2.  28073 
2.  60995 

2. 
2 

03586 
30707 
63257 

2.  08130 
2.  34581 
2.  66533 

2.  14598 
2.  40356 
2.71515 

2.  23029 
2.  48494 
2.  79069 

1  2.  33180 

2.  58999 

1  2.  89743 

2.  44775 

i  2.71620 
|  3.  U3485 

.9950 
.9975 
.9990 

2.  80703 

3.  02334 
1  29053 

2.  308S3 
Z  02500 

3.  29206 

2.  81432 

3.  03010 
3.  29673 

2.  83289 

3.  03S98 
3.  30489 

2.  S3830 

3.  05234 
3.31715 

I 
3- 

S5S94 
07144 
33464 

2.  S8859 

3.  09871 
3  35949 

2.  93347 

3.  13969 
3.  39647 

3.  00431 
3.  20556 
3.  45698 

1  3.  11073 

1  3.  31099 

3.  55939 

'  3.  25525 
1  3.  -»6164 
|  3.71692 

confidence  ellipse  is 


A  =  Ka  a  it 
e    x  y 


(2.38) 


where  K  is  the  appropriate  probability  conversion  factor 
(Table  I).   The  area  of  the  90  percent  confidence  circle  is 


Ac  =  ttR* 


(2.39) 


where  E  is  given  by  Equation  2.37.   For  a  condition  where 
o=a=3  meters,  and  8  =  30°,  the  area  of  the  90  percent 

1       2 

confidence  ellipse  is  261  square  meters,  while  the  area  or 
the  ccnfidence  circle  is  587  square  meters.   For  both  stan- 
dard errcrs  equaling  10  meters  and  8  =  30° ,  the  90  percent 
confidence  ellipse  has  an  area  of  921  square  meters  and  the 
confidence  circle  has  an  area  of  2894  square  meters.   From 
an  operational  perspective,  the  difference  in  areas  between 
ellipses  and  circles  have  significant  implications  which 
will  he  discussed  in  Chapter  V. 

The  following  examples  are  presented  to  demon- 
strate methods  for  computing  the  parameters  cf  error 
ellipses  ard  confidence  circles  for  several  hydrographic 
positioning  geometries. 
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Example  _1 

A  vessel  is  conducting  a  hydrographic  survey 

using  range-range  geometry.   The  two  LOP's  generated  ty 

microwave  transmitters  have  standard  errors  of  a  =  3  meters 

1 
and  a  =  4  meters.   The  angle  of  intersection  8  at  the 

2 

vessel    is   30°.       Assume    the    LOP's    are   uncorrelated.       Compute 
the   probability    that    the    vessel's    position    will    be    within   a 
circle    of    10-meter    radius   with    the    center    at    the    intersec- 
tion   of    the   LOP's. 

Eecalling   Equations   2.18    and    2.19,    the    values   of 

o      and    a     are    found    tc   be   9.79    meters   and    6.11   meters, 
x  y 

respectively.      From    Eguation    2.36 

a 

c   =    -*   =    0.633 
a 
x 

and  frcm  Eguation  2.37,  with  R  =  10  meters, 

K  =  1.032 
Entering  Table  III  and  using  interpolated  values  for  c  and 
K,  the  probability  that  the  vessel's  position  will  be  within 
a  circle  of  10-meter  radius  centered  at  the  intersection  of 
the  TCP's  is 

P  =  53.2% 


Example  2 

A  vessel  is  conducting  a  hydrographic  survey 
using  range-azimuth  geometry.   The  range  LOP  generated  hy 
the  microwave  transmitter  has  a  standard  error  of  3  meters. 
The  azimuth  LOP  determined  by  theodolite  observation  has  a 
standard  error  of  1.3  meters  at  all  ranges.   Compute  the 
radius  of  the  90  percent  confidence  circle  at  the  vessel*s 
positior. 

In  the  range-azimuth  case  p  =  90°  and  the  ICP's 

are  uncorrelated.   Therefore, 

a  =  a  =3.0  meters 
1      x 

and 
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a      =     a    =    1.  3    meters 


then  y 

a 
c   =   -i    =    0.433 
a 
x 

Table    IV   is   entered    with    the    values    of    ?    =    0.9    and    c   = 
0.433.      The   value   for   K   is    found    to    be 

K   =    1.7117 
Using    Equation    2.37,    the    radius    of    the    90    percent    prob- 
ability   circle    is    found    to   be 

R   =    5.14    meters 
The    probability    that    the    vessel's   position    will    be    within   a 
circle   of    5.14-meter    radius    centered    at    the    intersection    of 
the    ICP's   is   90    percent. 

Example    3 

A    vessel    is   conducting    a    hydrographic    survey 
using    hyperbolic-hyperbolic    geometry.       The    hyperbolic    LC? 
generated   by   the    1.6-MHz   electronic    positioning    system    has    a 
standarl    error    of    0.05-lane    on    the    base    line.       The    correla- 
tion   coefficient     (  p     )    between    the    two    LOP's    is    known    tc    be 

12 

0.4.       Compute   the    radius    of    the    90    percent    confidence   circle 
at    the    vessel's    position. 

The    rectangular    plane    coordinates  of    the    master 
(M)  ,    two   slaves     (G    and   P.),    and    the    vessel's   position     (?)    are 

I    COORDINATE  Y    COORDINATE 

(m)  (m) 

R  172,679.1  62,540.4 

G  308,679.1  98,540.4 

M  241,738.2  21,325.4 

P  223,172.5  169,264.2 
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Given    th€    frequency    of    1.6    MHz,    X    =    187.37    meters    frcn 

Equation    2.3.       The    lare    width    alcng    the    base    line    is  w1     =  w' 

g  r 

=   93.68    meters    from    Equation    2.5-      Using    the    law    of    cosines 

from    plare    geometry,    the    subtended   angles    a      and     a     are 

32.47°    and    43. 25°,    respectively.      The   angle    of   intersection 

of    the    tiio   hyperbolas  at    P   is    37.86°    from    Eguaticn    2.7.       The 

lane    widths   at    P   are    v      =    254.19    meters    and    w      =    335.06 

r  g 

meters    from   Equation    2.6.      The    standard    errors   of    the    green 

(a    )     and   red    (a  )    hyperbolas,    respectively  are    a  =    w  a           = 

12  i            g  base 

16.7    meters   and     a    =   w  a           =    12.7    meters.  These  standard 

2    r  base 

errors  are  in  a  linear  skewed  coordinate  system  and  must  be 

transformed  to  an  unccrrelated  rectangular  system.   Frcn 

Equations  2.18  and  2.19,  the  values  of  a  and  a      are  36.9 

a       b 

meters  and  12.7  meters,  respectively.   The  correlation  coef- 
ficient in  the  correlated  rectangular  system  (P  .  )  is  then 
0.737  from  Equation  2.26.   The  semi-major  and  semi-minor 
axes  in  the  uncorrelated  rectangular  system  are  38.1  meters 
and  8.3  meters,  respectively,  from  Eguations  2.27  and  2.26. 

The  eccentricity  is 

a 
c  =  ^  =  0.218 
x 

Table  IV  is  entered  liith  the  values  of  P  =  0.9  and  c  = 
0.218.   The  value  for  K  is  found  to  be 

K  =  1.6602 
From  Equation  2.37,  the  radius  of  the  90  percent  probability 
circle  is  found  to  be 

S  =  63.3  meters 
The  probability  that  the  vessel's  position  will  be  within  a 
circle  of  63.3-meter  radius  centered  at  the  intersection  of 
the  LCP's  is  90  percert. 
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Ill-  EXPERIMENT  DESIGN  AND  IMPLEMENTATION 

The  goals  of  this  chapter  are  to  demonstrate  that  hydro- 
graphic  positioning  accuracy  can  be  classified  based  or  the 
radii  cf  90  percent  confidence  circles  determined  by  using 
Eurt's  method  and  to  show  that,  based  on  the  same  criteria, 
accuracy  predictions  can  be  made  for  survey  planning 
purposes. 

A.   EATA  ACQUISITION  IBOCEDOBES 

The  data  used  for  analysis  and  prediction  consisted  of 
range-range,  azimuth-azimuth  and  range-azimuth  survey  infor- 
mation.  The  data  were  acquired  by  Naval  Postgraduate  School 
(NPS)  students  in  a  Hydrographic  Sciences  course.   Although 
the  ccurse  was  structured  as  a  training  exercise,  the  data 
acquisition  procedures  utilized  were  nearly  identical  to 
those  which  are  practiced  by  NOS. 

A  total  of  453  hydrographic  positions  were  recorded 
during  the  survey  of  a  nearshore  area  in  southern  Monterey 
Bay,  California.   Of  the  positions  used  for  analysis,  292 
were  range-range,  81  were  range-azimuth,  and  80  were 
azimuth-azimuth.   All  survey  information  was  recorded  by 
hand  in  sounding  volumes.   The  vessel  used  was  a  36-foot 
Uniflite  with  a  fiberglass  hull  and  twin  engines.   The 
survey  was  conducted  en  October  28,  November  16,  23,  and  30, 
1983.   Electronic  control  and  calibration  stations  used  for 
the  survey  included  CSE  MON  1978,  MOSSEL  1932,  BEACH  LAB 
1982,  MCNTEEEY  AMERICAN  CAN  COMPANY  STACK  1932,  MONTEREY 
RADIO  STATION  KMBY  MAST  1962,  MONTEREY  HARBOR  LIGHT  6  1978, 
and  MCNTEEEY  BLUE  LIGETHOUSE  (Fig.  3.1).   With  the  exception 
of  MCNTEEEY  BLUE  LIGHTHOUSE,  which  is  a  low-order  position, 
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MONTEREY   BAY 


121  '  55 


Figure  3.  1   Hydrographic  Survey  Area 

all  stations  are  of  third-order  or  better  and  are  published 
in  the  National  Geodetic  Survey  Data  Base. 

For  azimuth-azimuth  and  range-azimuth  positioning, 
azimuths  were  measured  with  a  Wild  T-2  theodolite.   Cn 
November  16,  range-azimuth  information  was  acquired  fcy 
locating  the  theodolite  over  station  MUSSEL  and  initialing 
cn  USE  MCN.   The  initial  direction  was  checked  by  sighting 
on  KMSY  MAST.   Azimuth-azimuth  positions  were  acquired  cn 
November  23.   A  theodclite  was  set  over  USE  WON  and  an 
initial  direction  was  to  MUSSEL   A  second  theodolite  was 
set  at  MOSSEL  using  USE  MON  for  the  initial  direction. 
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Range  information  was  recorded  using  a  Racal  Decca 
Trispcnder  systeir,  a  niicrowave  system  commonly  used  for 
nearshore,  line-of-sight  survey  work.   On  October  28  and 
Novemter  30,  range-range  data  were  recorded  by  setting 
remote  units  over  stations  BEACH  LAB  and  MUSSEL.   Eefore  and 
after  the  survey,  the  ranging  system  was  calibrated  ever  the 
fixed  tase  line  USE  MCN  to  MUSSEL.   Daily  checks  in  the 
survey  area  were  made  to  determine  if  the  system  was  working 
properly.   This  was  accomplished  by  maneuvering  the  survey 
vessel  tc  a  point  where  two  known  navigational  ranges  inter- 
sected.  One  navigational  range  was  formed  by  stations 
MONTEREY  AMERICAN  CAN  COMPANY  STACK  and  MONTEREY  RADIO 
STATICN  KMEY  MAST.   A  second  navigational  range  was  formed 
by  stations  MONTEREY  EABBOH  LIGHT  6  and  MONTEREY  BLUE 
LIGHTHOUSE. 

Track  control  for  range-azimuth  and  range-range  posi- 
tions was  accomplished  by  steering  the  vessel  along  range 
arcs.   The  spacing  between  range  arcs  for  most  lines  was 
planned  to  be  40  meters.   Distance  between  positions  along  a 
sounding  line  averaged  approximately  200  meters.   The 
azimuth-azimuth  lines  were  controlled  by  steering  a  magnetic 
compass  heading. 

The  data  acquired  under  training  conditions  contained 
several  deficiencies  that  would  normally  not  be  tolerated. 
For  example,  the  quality  of  the  line  steering  was  generally 
poor;  the  vessel  wandered  off  the  arc  more  than  10  meters  in 
several  instances.   The  quality  of  the  sounding  lines  run 
using  azimuth-azimuth  control  was  extremely  deficient;  the 
position  plot  of  these  lines  show  a  jagged  path  by  the 
vessel.   Under  normal  hydrographic  procedures,  these  posi- 
tions would  be  rejected.   Since  the  intent  of  this  study  is 
to  demonstrate  accuracy  analysis  techniques,  these  deficien- 
cies prove  to  be  inccrsequential;  the  acquired  data  are 
adequate  to  demonstrate  the  concepts. 
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IV.  RESULTS  AND  DATA  Mi^ISIS 

fi.   EATA  PROCESSING 

Automated  processing  of  the  positional  survey  data  was 
done  on  the  NPS  IBM  270/3033AP  computer  system.   Graphic 
displays  were  constructed  using  the  Display  Integrated 
Software  System  and  Plotting  Language  (DISSPLA)  developed  by 
the  Integrated  Software  Systems  Corporation  (ISSCO) 
[Ref.  16].   All  computer  programs  involved  in  data 
processing  were  written  in  the  WATFIV  programming  language. 

Computations  were  made  in  an  X-Y  coordinate  system  based 
on  a  Modified  Transverse  Mercator  (MTU)  projection.   A  MTM 
projection  is  essentially  the  same  as  a  Universal  Transverse 
Mercator  (UTM)  projection,  the  only  difference  being  that  in 
a  MTM  prcjection  a  central  meridian  is  picked  near  the 
survey  area  instead  of  being  fixed  at  a  particular  meridian 
[Ref.  17]. 

The  central  meridian,  controlling  latitude,  and  false 
easting  values  define  the  coordinate  system  used  for  compu- 
tations.  The  central  meridian  for  the  projection  was  chosen 
to  be  longitude  121°  52*  30"  W  which  is  approximately  the 
mean  longitude  of  the  survey  area.   The  controlling  lati- 
tude, tie  distance  iE  meters  from  the  eguator  to  a  reference 
latitude,  was  chosen  to  be  4,050,000  meters.   A  false 
easting  cf  5,000  meters  was  chosen  as  the  value  of  the 
X-coordinate  at  the  central  meridian. 

Three  shore  contrcl  stations  were  used  in  the  acquisi- 
tion of  survey  data.   The  geodetic  positions  of  these 
stations  were  converted  to  the  X-Y  coordinate  system  (Table 
V)  using  program  UCOMPS,  which  is  a  hydrographic  utility 
package  available  to  students  at  NPS. 
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TABLE  V 
Coordinates  of  Control  Stations 


STATION  NAMS 
DSE  HON 

rOSSEL 

EEACH  LAB 


GICDETIC  CCORD. 


360  36' 

1210  52' 

360  37' 

1210  5u» 

360  36' 

1210  52' 


04.685"  N 
35.900"  W 

18.151"  N 
11.628"  K 

05.571"  N 

33.427"  W 


MTfl  COORD. 

Y  =  1932.43  m, 
X  =  4853.36  it, 

Y  =  4247.42  m, 
X  =  2474.75  m, 

Y  =  2009.86  m. 
X  =  4914.75  m 


B.   ACCURACY  ANALYSIS  OF  HYDROGRAPHIC  POSITIONING  DATA 

The  objective  of  this  section  is  to  illustrate  how  the 
accuracy  of  hydrographic  positioning  data  can  be  classified 
using  Eurt's  method  of  circles  of  eguivalent  probability. 
The  radius  of  the  90  percent  confidence  circle  was  computed 
for  each  position;  it  provides  a  quantitative  measure  of 
repeatatle  accuracy. 

Fox  subsequent  accuracy  computations,  the  following 
assumptions  were  made: 

i.   The  standard  error  for  the  microwave  ranging 
system  used  in  the  range-range  and  range- 
azimuth  computations  is  3  meters. 
ii.   Eor  azimuth-azimuth  and  range-azimuth  positions, 
the  pointing  error  of  the  theodolite  is  1.3 
ireters  at  all  ranges, 
iii.   The  two  LOP's  involved  in  all  types 

cf  positioning  are  independent  (  p   =  0)  . 

1  2 

iv.   The  data  are  free  of  systematic  errors. 
Raw  range  and  azimuth  data  were  hand  logged  into  a  data 
file  for  processing.   A  modification  of  program  UCOMPS  was 
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used  to  compute  X-Y  coordinates  of  all  positions.  Based  en 
geometric  relationships  discussed  earlier,  angles  of  inter- 
section of  the  LCP's  viere  then  computed  for  range-range  and 
azimuth-azimuth  points.  The  angles  of  intersection  for  all 
range-azimuth  positions  are  90°. 

The  range-range  and  azimuth-azimuth  data  were  then 
passed  to  WATFIV  subroutine  PECB  (Appendix  A).   As  input 
parameters,  the  subroutine  accepts  two  standard  errors  of 
the  LCP's  and  the  corresponding  angle  of  intersection.   The 
output  parameters  include  the  semi-major  and  semi-minor  axes 
of  the  90  percent  confidence  ellipse,  the  radius  of  the  90 
percent  confidence  circle,  and  the  areas  covered  by  both 
figures. 

Subroutine  PEOB  uses  a  linear  approximation  to  determine 
the  value  of  the  function  K  for  varying  values  of  the  eccen- 
tricity, c,  in  Burt's  method.   A  linear  interpolation  was 
performed  by  first  taking  the  eleven  discrete  values  of  c 
and  K  for  a  probability  of  90  percent  from  Table  IV  and  then 
constructing  a  series  of  relationships  for  K  as  a  function 
of  c  (Table  VI) . 

Values  of  the  radii  of  90  percent  confidence  circles  for 
range-range  data  were  plotted  at  their  respective  positions 
(Fig.  4.1).   The  arcs  of  circles  connecting  the  two  control 
stations  BEACH  LAB  and  MUSSEL  represent  lines  of  constant 
intersection  angle  (30°).   Of  the  range-range  data  set, 
position  848  (Appendix  B) — coordinates  X  =  4119.01,  Y  = 
4735.  C7 — was  found  to  have  the  smallest  radius  (strongest 
position)  of  6.4  meters  and  an  angle  of  intersection  of 
90.20.   Position  137 — coordinates  X  =  3345.86,  Y  = 
3873-34 — represents  the  weakest  position  with  radius  value 
of  15.3  neters  and  an  angle  of  intersection  of  26.7°. 

The  positional  accuracy  degrades  rapidly  as  the  inter- 
section angle  approaches  30°;   the  30°  arc  represents  a  line 
of  constant  13.7  meter  radius.   Within  400  meters  of  the  30° 
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TABLE  VI 

Linear  Approxisations  for  K  as  a  Function  of 

c 

Interval  of  c 

Linear  Interpolation 
"Function  for  "K 

0.0  -  0.  1 

K  =  .0306c  +  1  .64485 

0.1  -  0.2 

K  =  .0940c  +  1.63851 

0.2  -  0.3 

K  =  .  1652c  +  1.62427 

0.3  -  0.4 

K  =  .2535c  +  1.59778 

0.4  -  0.5 

K  =  .3790c  ♦  1.54758 

0.5  -  0.6 

K  =  .5444c  +  1  .46433 

0.6  -  0.7 

K  =  .7101c  +  1.36546 

0.7  -  0.8 

K  =  .8508c  +  1.26697 

0.3  -  0.9 

K  =  . 9475c  +  1  .18961 

C.9  -  1.0 

K  =1. 0361c  +  1.10987 

intersection  arc,  the  radius  varies  between  8  and  15  meters. 
The  radii  values  charge  slowly  in  the  vicinity  of  the 
minimun  value  of  6.4  meters  which  corresponds  to  an  angle  of 
intersection  of  90°. 

The  radii  of  90  percent  confidence  circles  associated 
with  the  azimut h-azinuth  positions  acquired  using  control 
stations  USE  MON  and  F.USSEL  were  also  plotted-  at  their 
respective  positions  (Fig.  4.2).   The  standard  errors  of  the 
LOP's  are  assumed  to  te  1.3  meters;  the  resulting  improved 
accuracy  is  evident.   The  maximum  value  of  the  90  percent 
confidence  circle  radii  is  8.7  meters  at  position 
637 — coordinates  X  =  4327.25,  Y  =  2818.39 — which  corresponds 
to  an  angle  of  intersection  of  159.8°  (or  in  terms  of  the 
supplement,  20.2°).   Eosition  682 — coordinates  X  =  4611.20, 
Y   =  4421.29 — represents  the  strongest  position  recorded 
during  the  survey  with  a  90  percent  confidence  circle  radius 
of  2.8  meters  and  an  angle  of  intersection  of  91.0°. 

Again,  the  rapid  degradation  of  accuracy  is  noted 
approaching  8  =  150°.   The  arc  of  the  150°  intersection 
angle  represents  a  corstant  radius  of  5.9  meters.   Discrete 
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SURVEY  DATA  ANALYSIS  R/R 

RADII  OF  90%  PROBABILITY  CIRCLES 
STATIONS  BEACH  LAB  AND  MUSSEL 


3400 

EASTINGS 


Figure  4. 1    Raage-Hange  Accuracy  Analysis 
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Fiqure  4.2    A2imuth- Azimuth  Accuracy  Analysis 
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values  along  the  arc  confirm  this  qualitatively.   A  large 
area  of  strong  positional  accuracy  surrounds  the  area  where 
8   =  90°.   Numerous  values  of  2.8  meters  are  present  near 
the  top  cf  the  plot. 

Using  the  assumptions  stated  at  the  beginning  of  this 
section,  the  values  for  all  radii  of  90  percent  confidence 
circles  for  range-aziiuth  positions  are  5.1  meters.   This 
computation  was  carried  out  in  Example  2  of  Chapter  II. 
Since  this  case  is  trivial,  the  data  are  not  displayed 
graphically. 

Positioning  data  were  also  classified  based  on  the 
parameters  of  the  90  percent  confidence  ellipse.   WAIFIV 
program  F1LIP  (Appendix  C)  was  used  to  generate  the  parame- 
ters cf  the  90  percent  confidence  ellipse  for  range-range, 
azimuth-azimuth,  and  range-azimuth  positioning  data.   Ihe 
program  was  initialized  by  entering  the  coordinates  of  the 
control  stations  and  standard  errors  of  the  LOP's.   The  fix 
number,  hydrographic  position  coordinates,  and  angle  of 
intersection  were  then  read  in  from  a  data  file.   Subroutine 

PROB  vas   called  to  ccnpute  values  for  Kg   and  K a  . 

x         y 

The  angle  of  orientation  of  the  major  axis  of  the 
ellipse,  measured  clockwise  from  north,  was  then  computed. 
lor  range-range  and  azimuth-azimuth  positions,  the  LCF 
generated  from  the  left  control  station  was  used  as  the  base 
LOP.   For  range-azimuth  positions,  the  LOP  formed  by  the 
theodolite  was  used  as  the  base  LOP.   First,  the  orientation 
of  the  base  LOP  in  the  coordinate  system  was  determined. 
The  orientation  of  the  major  axis  of  the  error  ellipse  rela- 
tive to  the  base  LOP  (9)  was  then  computed  using  Eguaticn 
2.29.   By  adding  or  subtracting  9  to  the  orientation  of  the 
base  LOP,  the  orientation  of  the  major  axis  of  the  error 
ellipse  in  the  coordinate  system  was  determined.   This  angle 
takes  on  values  from  0°  to  180°.   Appendix  D  consists  cf  the 
confidence  ellipse  classification  scheme  for  range-range, 
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azimuth-azimuth  and  range-azimuth  data.   Forty  positions  for 
each  positioning  geometry  are  listed  for  comparison  to  the 
classification  scheme  presented  in  Appendix  3. 

Appendix  B  lists  the  data  by  position  number,  X-Y  coor- 
dinate, angle  of  intersection,  and  radius  of  the  90  percent 
confidence  circle.   Appendix  D  lists  the  data  by  position 

number,  X-Y  coordinate,  angle  cf  intersection,  Kg  ,  Kg  ,  and 

x     y 

angle  of  orientation  for  the  90  percent  confidence  ellipse. 
These  appendices  are  similar  to  hydrographic  survey  data 
bases  and  demonstrate  accuracy  classification  schemes  based 
on  the  two  criteria. 

C.   ACCURACY  PEEEICT1CHS 

lie  overall  positional  accuracy  of  a  survey  can  re 
controlled  by  computing  accuracy  values  before  data  acquisi- 
tion is  begun.   For  example,  if  the  hydrographer  is  using 
radii  cf  90  percent  confidence  circles  as  an  accuracy 
criterion,  the  minimum  allowable  angle  of  intersection  for 
two  LCP's  can  be  computed  for  meeting  specifications.   The 
nature  of  the  survey  area  may  allow  the  flexibility  to 
change  system  geometry  to  maximize  accuracy  a*t  a  specific 
location  or  to  maximize  the  area  covered  with  a  given  accu- 
racy.  By  making  accuracy  computations  before  acquiring 
data,  the  hydrographer  may  also  have  the  option  of  deciding 
what  type  cf  positioring  system  is  to  be  used  to  meet  accu- 
racy requirements. 

The  construction  cf  reliability  contours  is  one  method 
to  display  the  expected  positional  accuracy.   Reliability 
contours,  lines  cf  constant  repeatable  accuracy  which  are 
functions  of  the  system  geometry  and  standard  errors  cf  tie 
positicrirg  equipment,  can  be  constructed  about  shore 
stations  using  the  radii  of  90  percent  confidence  circles 
criterion  or  the  less  desirable  d^0  value. 
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Consider  the  equations  that  have  been  developed  in 
Chapter  II  for  the  determination  of  radii  of  90  percent 
confidence  circles  using  Burt's  method.   For  uncorrelated 
lOP's  in  a  range-range  or  azimuth-azimuth  system,  the 
repeatable  accuracy  cf  a  hydrographic  position  is  a  function 
only  cf  the  angle  of  intersection,  assuming  the  standard 
errors  of  the  LOP's  are  constant  throughout  the  survey  area. 
The  lccus  of  points  which  define  a  constant  angle  of  inter- 
section for  two  LOP's  in  a  range-range  or  azimuth-azimuth 
system  is  a  circle  which  passes  through  both  control 
stations.   Given  the  coordinates  of  the  two  control 
stations,  the  equations  of  these  circles  can  be  determined. 

Construction  of  reliability  contours  involves  several 
simple  trigonometric  relationships  (Fig.  4.3).   Let  IE  be 
the  line  connecting  the  two  shore  control  stations  L  and  E 
in  a  range-range  system.   The  length  of  line  LR  is  b.   The 
circle  through  both  stations  defines  a  line  of  constant 
intersection  angle  for  two  LOP's.   The  radius  of  the  circle 
is  r.   The  distance  e  is  measured  along  the  perpendicular 
bisector  of  the  line  IE  to  the  center  of  the  circle  at 
point  0(h,k)  and  is  given  by 


e  ■ 


2tan8  (4.1) 


Knowing  e  and  the  radius  r,  the  coordinates  of  point  0  can 
be  computed.   The.  equation  of  the  circle  is  then 


r1   «  (x  -  h)2  +  (y  -  k) 


(4.2) 


These  two  equations  were  used  to  generate  reliability 
contours  for  display  on  a  computer  graphics  terminal.   Using 
Eurt's  method,  the  angles  of  intersection  of  two  LOF's  were 
computed  for  discrete  values  of  radii  of  90  percent  confi- 
dence circles.   Reliability  contours  about  stations  EEACH 
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Figure  4.3    Ccistruction  of  a  Reliability  Curve 

LAB  and  HDSSEL  for  a  range-range  system  (a  =a  =  3  meters) 

1  2 

were    constructed    (Fig-    4.4).       Dsing    Equation    4.2,    X-Y    coor- 
dinates   Mere  generated   for  points   laying   on    different    reli- 
ability  circles.      A   curve- fitting   subroutine  in    the    EISSPLA 
library    was   used   to    generate   the   circles    through    the 
computed   points.      The    13-meter   accuracy   contour    corresponds 
to   an    angle   of    intersection   of   31. 6°,    while   the    7-ireter 
accuracy   contour  corresponds    to  an  angle   of   intersection   of 
67.9°.      The   best   achievable    accuracy    of    the    system    is    6.4 
meters    at   90°. 

For   comparison    purposes,    reliability  contours    were 
constructed   about    BEACH   LAB    and    MUSSEL   for    azimuth-azimuth 

geometry    (  a    =q     =1.3    meters) .      The    increased   accuracy    of 

i  2 

this    configuration    is    evident     (Fig.    4.5) .      The   3-meter 
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contour  corresponds  tc  an  angle  of  intersection  of  69.4° 
while  the  6-meter  contour  corresponds  to  an  angle  cf  inter- 
section of  29.6°.   The  best  achievable  accuracy  at  an  inter- 
section angle  of  90°  is  2.8  meters. 

A  second  scheme  *as  used  tc  display  accuracy  predictions 
for  the  twc  positioning  methods.   Given  the  coordinates  cf 
ESACH  LAB  and  MUSSEL,  a  series  of  discrete  points  spaced  800 
meters  apart,  were  generated  throughout  the  survey  area. 
The  values  for  the  radii  of  90  percent  confidence  circles 
were  then  computed  at  each  point  with  the  use  of  subroutine 
FROB.   Figures  4.6  and  4.7  illustrate  this  prediction 
scheme.   These  figures  present  the  same  information  as 
Figures  4.4  and  4.5  in  a  different  manner.   The  30°  angle  of 
intersection  contour  is  shown  on  both  figures. 
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GENERATED  RELIABILITY  CONTOURS 

RANGE-RANGE:  BEACH  LAB-MUSSEL 
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Figure    4.4        Reliability   Contours:       Range-Range    Geometry 
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GENERATED  RELIABILITY  CONTOURS 

AZIMUTH-AZIMUTH:  BEACH  LAB-MUSSEL 
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Figure  4.5   Beliability  Contours:   Azimuth- Azimuth  Geometry 
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ACCURACY  PREDICTIONS  RA-RA 

RADII  OF  90%  PROBABILITY  CIRCLES 
STATIONS  BEACH  LAB  AND  MUSSEL 
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Jigure  4.6   Hange-Hange  Point  Accuracy  Prediction 
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ACCURACY  PREDICTIONS  AZ-AZ 

RADII  OF  90%  PROBABILITY  CIRCLES 
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Figure   4.7        Azimuth-Aziauth   Point   Accuracy   Prediction 
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V.  CONCIDSIONS  AND  RECOMMENDATIONS 

A.   ACCURACY  SPECIFICATIONS 

Interpretation  of  the  1982  IHO  positioning  standards  in 
terms  cf  90  percent  confidence  circles  yields  some  inter- 
esting results  with  respect  to  present  day  survey  practices. 
For  example,  for  a  1  :  10,000-scale  hydrographic  survey,  NOS 
usually  uses  microwave  positioning  systems  in  a  range-range 
mode,  and  assumes  a  standard  error  of  3  meters  for  each  LCP. 
Surveys  are  freguently  conducted  between  the  30°  to  150° 
angle  of  intersection  limits.   Osing  the  90  percent  confi- 
dence circle  criterion,  the  radius  of  the  circle  should  net 
exceed  10  meters.   However,  the  radius  value  for  8  =  20°    and 

150°  is  13.7  meters.   The  values  of  K a  and  Kg   for  the  90 

x       y 

percent  confidence  ellipse  are  17.6  and  4.7  meters,  respec- 
tively.  To  meet  the  SO  percent  criterion  for  a  1:  10,000- 
scale  survey,  the  6  limits  should  be  42°  to  138°. 

Azimuth-azimuth  positioning  is  accurate  enough  for 
1: 5 ,0 OC-scale  surveys,  using  8  limits  of  35°  to  145°, 
assuming  a  standard  error  of  1.3  meters  for  each  LCP.   Kith 
the  standard  error  assumptions  used  for  range-azimuth,  the 
90  percent  radius  is  5.1  meters  for  all  positions.   Given 
the  uncertainties  of  the  standard  error  figures,  it  is 
rational  to  assume  that  range-azimuth  positions  can  meet  the 
5-meter  accuracy  standard  for  1: 5, 000-scale  surveys.   In 
fact,  range-azimuth  positional  accuracy  can  exceed  azimuth- 
azimuth  accuracy  when  the  later's  8  is  less  than  35°.   For  a 
3-meter  a  range-range  configuration,  it  is  impossible  to 
meet  1:5,000  specifications  with  any  8. 

As  a  general  guideline,  the  30°  to  150°  angle  of  inter- 
section limit  is  a  good  rule  to  use  for  uncorrelated  LCE's. 
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However,  as  mentioned  for  1 : 10 ,000-scale  surveys  in  a  range- 
range  mode  (a  =  3  meters),  this  rule  does  not  always  hold. 
On  the  other  hand,  it  is  possible  to  have  6  's  of  less  than 
30°  and  still  meet  specif icaticns.   For  example,  azimuth- 
azimuth  positioning  can  theoretically  be  used  for  8's  of  18° 
to  162°  for  a  1 :  10, OOC-scale  survey.   However,  the  eccen- 
tricity cf  the  error  ellipse  is  so  small  that  the  distortion 
introduced  by  using  confidence  circles  can  become 
misleading.   In  view  cf  this,  eccentricities  of  less  that 
0.2  should  not  be  used. 

Using  the  90  percent  radius  criterion,  a  table  has  teen 
assembled  illustrating  the  8  limit  for  various  positioning 
geometries  at  different  survey  scales,  using  assumed  stan- 
dard errcrs  (Table  VII) .   The  information  in  Table  VII 
illustrates  that  the  20°  to  150°  8  limit  need  not  be  fixed. 
The  8  limits  should  vary  based  on  the  scale  of  the  survey 
and  the  precision  of  the  positioning  eguipment.   Accuracy 


Survey 
Scale 


1:2,500 

1:5,000 

1:  IC'OOO 

1:2C,000 

1  :40,000 


TABLE  711 
3  Limits  for  Surveys 


90% 
Radius 


(m) 

2.5 

5.0 
10.0 
20.0 
40.0 


(deg) 


42-138 
27-153 
23-157* 


R-R 

a  =  10) 

Limit 

(deg) 


3  5-145 


Az-Az 
(  a  =  1.3) 
8  Limit 

(deg) 


35-145 
23-157* 
23-157* 
23-157* 


*  Eccentricity  limit  of  0.2 

Note:   90%  radii  cf  all  range-azimuth  positions  are 

assumed  to  be  5.1  meters  for   a  =  3   and  a   =  1.3. 

1  2 
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figures    as    a   function  of   8    for    uncorrelated   LCP's    have    teen 
compiled   using    standard   errors   of    1.3    meters   for    azimuth- 
azimuth     (Table    VIII)  ,3.0    meters    for    range-range    short-range 
(Table   IX) ,    and    10    meters    for    range-range    medium-range 
(Table    X)     positioning  systems. 


TABLE  7III 

Accuracy  Figures  for  a     =  a  = 

1.3  m,   d   =0 

1       2 

12 

Angle  of   Ka 
Inter.       x 

K  a     Radius  of 

Area  of    Area  of 

y     90%  Circ. 

Ellipse    Circle 

(deg)      (m) 

(m)       (m) 

(sg  m)      (sq  m) 

90 

2.8 

2.8 

2.8 

24    1     24 

85 

2.9 

2.7 

2.8 

25         25 

80 

3.  1 

2.6 

2.8 

25    |     25 

75 

3.2 

2.5 

2.9 

25    |     26 

70 

3.4 

2.4 

3.0 

26    |     28 

65 

3.7 

2.4 

3.  1 

27    |     30 

60 

3.  9 

2.3 

3.3 

28    |     34 

55 

4.7 

2.3 

3.5 

30    |     38 

50 

5.  1 

2.2 

3.8    j 

32    |     44 

45  . 

5.  1 

2.  1 

4.  1 

35    |     53 

40 

5.8 

2.  1 

4.5 

.38    |     65 

35 

6.6 

2.  1 

5.  1 

43         83 

30 

7.6 

2.0 

5.9 

49        110 

25 

9.  1 

2.0 

7.1 

58    |    156 

20 

11.4 

2.0 

8.8 

71    |    241 

15 

15.  1 

2.0 

11.6 

94        425 

10 

22.6 

2.0 

17.4 

141    I    949 
281      3,781 

5 

45.2 

2.0 

34.7 

K  =  2.146  for  ( 

30^  probability 
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...      — 

TABLE    II 

accuracy   Figores   for    a     =  o 

12 

1   3    a,  o 

12 

=    0 

Angle   of          Ka              Ko           Radius   of 

Area   of 

Area    cf 

Inter. 

*                7         90S   Circ. 

Ellipse 

Circle 

(deg) 

(n)              (n)                 (m) 

(sg    a) 

(sg    in) 

90 

6.4 

6.4 

6.4 

130 

1            130 

85 

6.7 

6.2 

6.5 

131 

131 

80 

7.  1 

6.0 

6.5 

132 

135 

75 

7.5 

5.7 

6.7 

135 

141 

70 

7.9 

5.6 

6.9 

139 

149 

65 

8.5 

5.4 

7.2 

144 

162 

60 

9.  1 

5.2 

7.5 

150 

179 

55 

9.9 

5.  1 

8.0 

159 

203 

50 

10.8 

5.0 

8.7 

170 

235 

45 

11.9 

4.9 

9.4 

184 

280 

40 

13.3 

4.8 

10.5 

203 

345 

35 

15.  1 

4.8 

11.8 

227 

440 

30 

17.6 

4.7 

13.7 

260 

588 

25 

21.0 

4.7 

16.3 

308 

832 

20 

26.2 

4.6 

20.2 

381 

1,284 

15 

34.9 

4.6 

26.8 

503 

2,262 

10 

52.2 

4.6 

40.  1 

750 

5,052 

5 

104.  4 

4.6 

80.  1 

1,494 

20,  135 

K   =    2.146    for    90S   probability 

TABLE    I 

Accuracy    Fi 

gures   for    o     =  o 

1              2 

=      10      U,         p 

12 

=    0 

Angle   of        Ko 

Ko            Radius    of 

Area   of 

Area   of 

Inter. 

X 

y          90S    Circ. 

Ellipse 

Circle 

(deg) 

(m) 

(a)                  (m) 

(sg    a) 

(sg    a) 

90 

21.5 

21.5 

21.5 

1,447 

1,447 

85 

22.5 

20.6 

21.6 

1,  452 

1,455 

80 

23.6 

19.8 

.   21.8 

1,469 

1,497 

75 

24.9 

19.  1 

22.3 

1,  498 

1,562 

70 

24.5 

18.5 

23.0 

1,540 

1,657 

65 

28.2 

18.0 

23.9 

1,596 

1,798 

60 

30.3 

17.5 

25.2 

1,671 

1,989 

55 

32.9 

17.1 

26.8 

1,  766 

2,252 

50 

35.9 

16.7 

28.9 

1,889 

2,615 

45 

39.6 

16.4 

31.5 

2,  046 

3,117 

40 

44.4 

16.1 

34.9 

2,251 

3,835 

35 

50.5 

15.9 

39.5 

2,522 

4,890 

30 

58.6 

15.7 

45.6 

2,894 

6,528 

25 

70.1 

15.5 

54.3 

3,  423 

9,250 

20 

87.4 

15.4 

67.4 

4,230 

14,270 

15 

116.3 

15.3 

89.4 

5^590 

25,126 

10 

174.  1 

15.2 

133.7 

8,332 

56,130 

5 

347.9 

15.2 

266.9 

16,600 

223,727 

k  =  ; 

2.146    for    ! 

)0J    probability 
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B.   DSES  FOB  ACCURACY  FIGURES 

NCS  is  currently  developing  the  Shipboard  Data  System 
III  (SDS  III) ,  a  hydrcgraphic  data  acquisition  and 
processing  system  which  will  replace  the  present  HYDROLCG/ 
HYDRCfLOI  system.   SCS  III  will  revolutionize  data  acquisi- 
tion and  processing  techniques  with  the  capability  tc 
perforn  high-speed  calculations  and  display  color  graphics. 
With  this  increased  computer  potential,  data  manipulat icns-- 
such  as  accuracy  computations—can  be  performed. 

Each  position  in  a  survey  can  be  given  a  quality  figure 
based  on  the  radius  of  the  90  percent  confidence  circle. 
This  figure  is  sufficient  for  non-critical  positions  of 
ordinary  hydrographic  data.   Critical  positions  are  these 
which  are  determined  for  significant  features  (i.e.,  wrecks, 
least  depths,  rocks,  and  other  potential  hazards).   Fcr 
these  positions,  the  parameters  of  the  90  percent  error 
ellipse  can  be  computed,  as  well  as  the  radius  of  the  90 
percent" confidence  circle. 

Many  schemes  can  le  envisioned  for  the  use  of  an  accu- 
racy figure-   For  exanple,  suppose  the  position  of  a 
submerged  pile  was  determined  by  range-azimuth  geometry  in  a 
prior  survey.   The  radius  of  the  90  percent  confidence 
circle  is  then  5.1  meters  (Ex.  2,  Ch.  II).   The  charting 
agency  now  wishes  to  relocate  the  pile  to  determine  if  it 
still  exists  and  is  still  a  hazard  to  navigation.   In  lew 
water  visibility,  a  cemmon  technique  used  to  resolve  such  an 
item  hould  be  to  send  divers  down  over  the  reported  position 
and  conduct  a  circle  search.   One  diver  remains  at  the 
reported  position,  holding  a  line,  while  the  other  diver 
swims  a  circumference  holding  the  other  end  of  the  line. 
Theoretically,  if  the  line  is  about  5  meters  long  and  a  hang 
does  net  occur,  it  is  90  percent  certain  that  the  pile  has 
been  removed.   For  a  higher  confidence,  the  line  is 
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lengthened.   In  an  investigation  such  as  this,  it  is  advi- 
sable tc  he  conservative  and  use  the  maximum  length  of  line 
which  is  operationally  feasible  to  provide  coverage  cf  an 
area  as    large  as  possible.   The  radius  of  the  90  percent 
confidence  circle  gives  the  hydrographer  a  rough  figure  fcr 
answering  the  guesticn:   Does  the  submerged  pile  exist? 

Knowing  the  parameters  of  the  error  ellipse  could  te 
useful  fcr  conducting  wire-drag,  wire-sweep,  and  side  scan 
sonar  operations.   Fcr  a  position  obtained  with  low  preci- 
sion positioning  eguipment,  the  search  to  relocate  a 
submerged  feature  could  cover  a  large  area.   Knowing  the 
parameters  cf  the  error  ellipse  could  reduce  the  area,  time, 
and  effort  of  the  search.   The  search  pattern  could  be 
planned  to  cover  the  desired  confidence  ellipse. 

With  the  guantif ication  of  accuracy,  a  decision  must  be 
made  concerning  how  much  confidence  is  needed  to  delete  a 
certain  feature  from  the  chart  after  a  search  has  been  made. 
The  90  percent  confidence  level  may  be  too  low,  whereas  the 
S5  or  99  percent  level  may  suffice.   A  balance  must  te  main- 
tained between  confidence  of  disproval  and  time  and  effort 
spent  on  the  search. 

Accuracy  predictions  in  the  form  of  reliability  contours 
can  be  displayed  using  computer  graphic  terminals.   These 
displays  will  contribute  to  the  efficient  planning  of 
surveys  to  meet  specifications.   Given  the  survey  area,  the 
available  ccntrol,  the  positioning  methods,  and  the  preci- 
sion of  the  positionirg  eguipment,  the  hydrographer  can  plan 
the  accuracy  of  the  survey  before  it  is  conducted.   The 
survey  area  and  the  available  control  may  be  such  that  there 
is  flexibility  to  change  control  stations  to  optimize  accu- 
racy over  an  area  of  critical  importance.   This  information 
can  be  displayed  graphically  and  plans  for  the  survey  can  be 
made  accordingly.   Likewise,  given  an  accuracy  limit,  such 
as  a  10-meter  radius  cf  the  90  percent  confidence  circle, 
the  area  tc  be  covered  at  that  accuracy  can  be  maximized. 
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Many  variables  exist  when  considering  accuracy  require- 
ments fcr  a  hydro-graphic  survey.   In  general,  higher  accu- 
racy means  more  time,  money,  and  effort.   Azimuth-azimuth 
geometry  is  the  most  accurate  method  of  positioning  analyzed 
in  this  thesis.   This  method  involves  at  least  two  people 
ashore  and  good  ship-to-shore  communications.   Currently, 
NOS  acquires  these  data  manually,  which  minimizes  the  speed 
that  the  vessel  can  operate  and  adds  to  processing  time.   Or. 
the  other  hand,  a  survey  using  a  medium-range  system  needs 
little  shcre  support  and  the  data  acquisition  is  automated. 
Accuracy  predictions  help  keep  a  balance  between  accuracy 
and  effort.   If  the  desired  accuracy  is  attainable  using  a 
range-range  system  instead  of  an  azimuth-azimuth  system, 
then  the  chcice  is  ohvious. 

Hydrcgraphic  positioning  in  the  future  will  be  dominated 
hy  two  methods.   For  cffshore  surveys,  the  Global 
Positioning  System  {GIS)    is  expected  to  give  positional 
accuracy  to  10  meters  or  better.   GPS  is  a  satellite  posi- 
tioning system  currently  being  deployed  by  the  Department  of 
Defense  and  will  provide  near  worldwide  coverage  for  users. 
Since  the  full  constellation  of  18  satellites  will  not' he 
operational  until  1988,  it  is  not  yet  known  if  the  expected 
accuracy  of  10  meters  will  be  met.   Nearshore  surveys  may 
use  multiple  LOP's  for  establishing  hydrographic  positions. 
The  principle  of  least  squares  is  applied  to  redundant 
observations  yielding  the  most  probable  position.   Fcr  roth 
GPS  and  least  squares  positioning,  confidence  ellipses  and 
circles  can  be  determined,  although  the  techniques  involved 
are  much  mere  complicated  than  those  presented  in  this 
thesis. 

The  accuracy  classification  scheme  presented  in  this 
thesis  is  predicated  en  the  elimination  of  systematic 
errors.   Much  work  is  needed  in  identifying  the  sources  of 
systematic  errors  associated  with  hydrographic  positioning 
equipment. 
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APPENDIX  A 
SUEBCUTINE  FOE  90  PERCENT  CONFIDENCE  CIRCLE  PARAMETEES 


SOEECUTINE  PEOE  (SIG 1 , S 1G 2 ,COR , TBETA, SGX90, SGY9  0, 

*  RADIUS, ELAR.CIEAR) 
IMPLICIT  REAL*4  (A-H.O-Z) 

C  COMPOTES  RADIOS  OF  90%  CONFIDENCE  CIRCLE  (BURT,  METHOD  2) 

C  THIS  SOBROUTINE  WORKS  FOE  COERELATED  AND  UNCOERELATED 

C  LINES  OF  POSITION. 

C 

C  INPUT  PABAMETEES: 

C         SIG1  AND  SIG2-  STANDARD  ERRORS  OF  TWC  LOP'S 

C         TBETA  -  ANGLE  OE  INTERSECTION  IN  DEGREES 

C  (0-180  DEG) 

C        COR  -  CORRELATION  COEFFICIENT 

C  (USUALLY  ZERO  EXCEPT  FOR  HYPEEECIIC 

C  OR  SEXTANT  POSITIONING) 

C  O0TE0T  PAEAMETERS: 

C        SGX90  AND  SGY90-  SEMI-MAJOR  AND  MINOR  AXES 

C  OF  90*  ERROR  ELLIPSE 

C        RADIUS  -  RADIOS  OF  90%  CONFIDENCE  CIRCLE 

C        ELAR  -  AREA  CF  90%  CONFIDENCE  ELLIPSE 

C        CIRAR  -  AREA  OF  90%  CONFIDENCE  CIRCIE 

C 

C   WORK    WITH    AN    ANGLE    IESS    THAN    90    DEGREES 

IF  (TEETA    .GT.90.)     BETA    =180. -TBETA 

IF  (TBETA. LE. 90.)     BETA=TBETA 
C    CHANGE    DEGREES    TO    EADIANS 

RAE=.0174532*BEIA 
C    TRANSFORMATION    SIG1    AND    SIG2    TO    CORRELATED      ' 
C    RECTANGULAR    SYSTEM 

SIGA=SQRT  (  (1-/<  (SIN  (RAD)  )  **2)  )  *  (SIG  1**2+    2.  *C0R*SIG1* 
*SIG2*COS  (RAD)  +SIG2**2) -SIG2**2) 

^TGR=t!Tr/ 
C    TRANSFORM      CORREIATICN    COEFFICIENT    TO    CORRELATED 
C    RECTANGULAR    COORDINATE    SYSTEM 

A=]  (SIG2*COS  (EAE)  J/SIG1)  +COR 

F=1/SQRT (1+2*CCE*SIG2*COS (RAD) /SIG1+  (SIG2/SIG1)  **2* 
* (CCS (RAD)  )  **2) 

COBAE=A*F 
C  TRANSFORM  TO  UNCORRELATED  RECTANGULAR 

AA=SQET ( (SIGA**2+SIGB**2)/2) 

CC=SQRT]l-  (4*SIGA**2*SIGB**2*  (1-CORAB**2) )/ 

*  (SIGA**2+SIGB**2)  **2) 
DD=SQRT (1+CC) 

SIGX=  AA*DD 

SIG Y  =  SQRT  (SIGA**2+SIGB**2-SIGX**2) 
C 
C  COMPUTE  ECCENTRICITY  OF  ELLIPSE 

C=SIGY/SIGX 
C  COMIUTE  BURT'S  K  FACTOR  BY  LINEAR  INTERPOLATION 
IF  ]  (CLE.  1.)  .AND.  (CGT.0.9)  )  THEN 
B=T.0361*C+ 1.10987 
ELSE  IF   ((CLE  .  0  .  9)  .  AND.  (C  GT.  0  .  8)  )  THEN 

B=0.9475*C+1. 18961 
ELSE  IF  (  (CIE-0.8)  -AND.  (CGT.  0.7)  )  THEN 

B=0. 8508*6*1.26697 
ELSE  IF  ((CIE.0.7)  .AND.  (C.GT.O.  6)  )  THEN 
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B=0. 7101*C+1.36546 
ELSE  IF  (  (C.IE.0.6)  .AND.  (C.GT.0.5)  )  THEN 
B=C. 5444*0+1.46488 

ELSE  IF  (  (C.IE.0.5)  .AND.  (C.GT.0.4) )  THEN 

B=0.3790*C+1. 54758 
ELSE  IF  {  (C.IE.0.4)  .AND.  (C.GT.0.3)  )  THEN 

B=0.2535*C+  1.59778 
EISE  IF  {  (C.IE-0.3)  .AND.  (C.GT.0.2) )  THEN 

B=0.1652*C+1. 62427 
ELSE  IF  (  (C.IE.0.2J..AND.  JC.GT.0.1))  THEN 

B—  .  r 
EISE 


,094*C+1. 63851 
B=0.0306*C+1. 64485 


END  IF 

EADIUS=B*SIGX 

SGX90=2. 146*SIGX 

SGY90=2. 146*SIGY 

CIRAR=3.  1415926*RADIUS**2 

EIAR=3.  1415926*SGX90*SGY90 
BEIUKN 
END 
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APPENDIX    B 
ACCURACY    CLASSIFICATION:       90    PERCENT    CONFIDENC1    CIRCUS 

CLASSIJIED    RANGE-RANGE    POSITIONS 
Control    Stations:      BEACH    LAB    1982    and   MUSSEL    1932 
Standard    Error    Used    in   Computations:    3    meters 


Fix 

X 

Y 

Angle  of 

Radius  of 

No. 

Coordinate 

Coordinate 

Intersection 

111    Circle 

1 

2668.05 

4942.07 

127.0 

8.2 

2 

2852.01 

5076. 69 

121.6 

7.7 

3 

3041.58 

5194.85 

113.6 

7.4 

4 

3040.27 

5148. 74 

117.0 

7.3 

5 

2838. 26 

5021.35 

120.2 

7.6 

6 

2640.06 

4877.  13 

126.9 

8.2 

7 

2553.46 

4752.36 

130.4 

8.7 

8 

2724.  10 

4885.47 

121.3 

7.7 

9 

2889.06 

5013.63 

117.6 

7.3 

10 

3075.61 

5124.37 

115.0 

7.2 

11 

3172.63 

5136.23 

112.7 

7.0 

12 

2958.  13 

5006.65 

114.  4 

7.  1 

13 

2771.33 

4876.02 

117.9 

7.4 

14 

2581. 02 

4729.89 

126.9 

8.2 

15 

2584.41 

4665.31 

124.0 

7.9 

16 

2740.05 

4805.63 

116.7 

7.3 

17 

2913.27 

4928.51 

112.8 

7.0 

18 

3097.44 

5047.77 

111.2 

6.9 

19 

3193.94 

5103.46 

110.9 

6.9 

20 

2627.47 

4659.05 

118.8 

7.4 

21 

2762. 15 

4769.60 

113.2 

7.  1 

22 

2904.91 

4875-24 

110.5 

6.9 

23 

3056.38 

4975.73 

109.3 

6.9 

24 

3207.36 

5064.00 

108.9 

6.8 

25 

3373.36 

5101.96 

107.  1 

6.8 

26 

3190. 64 

5007.38 

106.8 

6.8 

27 

3015.45 

4900. 75 

107.  1 

6.8 

28 

2839.38 

4778.31 

108.7 

6.8 

29 

2679.99 

4651.31 

112.8 

7.0 

30 

2727.88 

4632.77 

106.9 

6.8 

31 

2699.83 

4776.97 

105.2 

6.7 

32 

3092.77 

4901.75 

104.4 

6.6 

33 

3295.42 

5019.92 

105.0 

6.7 

34 

3502. 07 

5123. 36 

106.0 

6.7 

35 

3697.02 

5151.09 

105.3 

6.7 

36 

3474.77 

5065. 57 

104.0 

6.7 

37 

3257.48 

4953.37 

102.7 

6.6 

38 

3043.03 

4821.88 

101.7 

6.6 

39 

2845.07 

4680.80 

101.7 

6.6 

40 

2746.96 

4608.32 

103.  1 

6.6 

41 

2748. 09 

4550.40 

97.5 

6.5 

42 

2931.78 

4701.21 

98.4 

6.5 

43 

3134.33 

4834.56 

99.4 

6.5 

44 

3317.94 

4944.55 

101.0 

6.6 

45 

3515.  14 

5043.56 

102.6 

6.6 

46 

3720. 18 

5117.  15 

103.9 

6.6 

48 

3724. 93 

5083. 62 

102.6 

6.6 

49 

3514.33 

4989.  13 

100.3 

6.6 

50 

3312. 88 

4891.55 

98.5 

6.5 
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RANGE-RAKGE    ACCURACIES     (CONTINUED) 


Fix 

X 

Y 

Angle  of 

Radius  of 

No. 

Coordinate 

Coordinate 

Intersection 

907c  Circle 

51 

3114.26 

4776.58 

96.5 

6.5 

52 

2921.85 

4647.63 

94.8 

6.5 

53 

2753.76 

4502.52 

91.5 

6.4 

54 

2793.3  1 

4487.76 

86.6 

6.4 

55 

2960.81 

4619.38 

90.6 

6.4 

56 

3160.  14 

4750.25 

93.6 

6.4 

57 

3325.  11 

4847.99 

96.0 

6.5 

58 

3492. 32 

4942.05 

98.4 

6.5 

59 

3693.83 

5032.60 

100.8 

6.6 

61 

3568.42 

4910.65 

96.3 

6.5 

62 

3370.43 

4812.54 

93-4 

6.4 

62 

3174.72 

4695.41 

89.7 

6.4 

64 

2998.35 

4566.29 

84.5 

6.5 

65 

2824.36 

4431.  37 

76.9 

6.6 

66 

2819.  19 

4368.35 

67.7 

7.0 

67 

3002. 82 

4519.53 

80.0 

6.5 

68 

3184. 74 

4657.79 

86.  9 

6.4 

69 

3378. 02 

4763.00 

90.5 

6.4 

70 

3578.71 

4871.32 

94.4 

6.5 

71 

3788.85 

4962. 48 

97.7 

6.5 

72 

4010.25 

5031.39 

100.4 

6.6 

73 

3732.84 

4895. 17 

95.0 

6.5 

74 

3509.99 

4799.08 

91.3 

6.4 

75 

3307.  18 

4688.48 

86.9 

6.4 

76 

3094.50 

4546.  29 

80.  1 

6.5 

77 

2901.06 

4391.68 

68.5 

7.0 

78 

2754.33 

4263.59 

49.5 

8.7 

79 

3754. 15 

5020.93 

100.  1 

6.5 

80 

3514.53 

4919.46 

97.2 

6.5 

81 

3287.25 

4788.27 

93.3 

6.4 

82 

3062.52 

4651.  12 

89.4 

6.4 

84 

2846.46 

4482.96 

82.4 

6.5 

84 

2845.  13 

4562.17 

91.3 

6-4 

85 

2995.88 

4403-69 

68-0 

7.0 

86 

3165. 76 

4254.99 

52.7 

8.3 

87 

3331.61 

4102.46 

43.3 

9.8 

88 

3500.54 

3946.44 

37.5 

1  1.  1 

89 

3672.90 

3785.70 

34.0 

12.2 

90 

3846.42 

3623.  19 

32.0 

12.9 

91 

4022.93 

3464.  18 

31.6 

13.0 

92 

4406.76 

3998.00 

68.3 

7.0 

93 

4263.63 

4159.42 

70.3 

6.9 

94 

4112. 66 

4331.20 

73-9 

6.7 

95 

3959.57 

4527.77 

79.9 

6.5 

96 

3810.03 

4681.47 

85.5 

6.5 

97 

3672.32 

4864.  18 

93.7 

6.5 

98 

3535.90 

5040.98 

102.3 

6.6 

100 

3065.44 

4352.97 

61.8 

7.4 

101 

3068.28 

4166.32 

41.6 

10.1 

102 

3077.08 

4046.24 

29.5 

13.9 

103 

3117.61 

4262. 88 

52.8 

8.3 

104 

3062.  15 

4499.  1  1 

76.5 

6.6 

105 

3011. 82 

4658.04 

91.7 

6.4 

106 

3122.91 

4453.24 

71.4 

6.8 

107 

3152.  10 

4199.54 

47.  1 

9.  1 

108 

3118. 73 

4025.84 

29.3 

1  4.0 

110 

3188.79 

4154. 61 

43.8 

9.7 

11  1 

3177.20 

4396.30 

65.9 

7.  1 

112 

3111-79 

4583.07 

82.8 

6.5 

1  13 

3162.67 

4534. 87 

77.6 

6.6 

114 

3099.29 

4691.03 

91.3 

6.4 

115 

3209.60 

4484.20 

73.3 

6.7 

116 

3235.42 

4271.  54 

55.2 

8.0 

117 

22C4. 85 

4043.52 

34.3 

12.0 
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RANGE-RANGE    ACCURACIES      (CONTINUED) 


Fix 

X 

Y 

Angle  of 

Radius  of 

No. 

Coordinate 

Coordinate 

Intersection 

90S  Circl€ 

120 

3261.07 

4111.  17 

42.0 

10.0 

121 

3262.51 

4368.76 

63.7 

7.3 

122 

3175.94 

4613.63 

83.8 

6.5 

123 

3189.26 

4667.90 

87.5 

6.4 

124 

3286.57 

4471.03 

71.9 

6.8 

125 

3317.76 

4240. 01 

53.9 

8.  1 

126 

3282.62 

4020. 86 

35.3 

1  1.7 

127 

3251.48 

3951.93 

28.6 

14.3 

128 

3309.62 

3962.87 

31.8 

13.0 

129 

3362.57 

4200.72 

51.7 

8.4 

130 

3331.06 

4454. 74 

70.7 

6.9 

131 

3252.99 

4654.02 

85.4 

6.5 

132 

3286.69 

4673-74 

86.3 

6.5 

133 

3366.26 

4474.82 

72.2 

6.8 

134 

3387.38 

4275.01 

57.7 

7.8 

135 

3376.70 

4076.27 

42.6 

9.9 

136 

3337.33 

3910.52 

29.0 

1  4.  1 

137 

3345. 86 

3873.34 

26.7 

15.3 

138 

3422.73 

4077.78 

44.0 

9.6 

139 

3431.43 

4314. 80 

61.3 

7.4 

140 

3396.21 

4530„82 

76.0 

6.6 

141 

3326.96 

4700.26 

87.4 

6.4 

142 

3361.30 

4694„84 

86.7 

6.4 

143 

3436.64 

4506- 04 

74.4 

6.7 

144 

3474.89 

4309.23 

61.5 

7.4 

145 

3464. 73 

4120.80 

48.2 

8.9 

146 

3425.35 

3933. 68 

34.0 

12.2 

147 

3459.34 

3928.23 

34.8 

1  1.9 

148 

3513.95 

4137.04 

50.6 

8.5 

149 

3512.94 

4390.84 

67.4 

7.0 

150 

3473.67 

4567. 85 

78.4 

6.6 

151 

3385. 01 

4766.90 

90.7 

6.4 

152 

3390.09 

4801.65 

92.5 

6.4 

153 

3472. 78 

4657.67 

83.8 

6.5 

154 

3527.62 

4492.63 

73.9 

6.7 

155 

3555.04 

4377.14 

67.0 

7.  1 

156 

3561. 04 

4288.00 

61.4 

7.4 

157 

3543.95 

4204.99 

55.7 

7.9 

158 

3543.69 

4262.39 

59.5 

7.6 

159 

3543.85 

4068.09 

46.8 

9.  1 

160 

3493.62 

3880.  18 

32.9 

12.5 

719 

3519.28 

3843.16 

31.6 

13.0 

720 

3586.31 

4061.28 

47.6 

9-0 

721 

3589.75 

4237.59 

58.7 

7-7 

722 

3568.78 

4441.  15 

71.  1 

6.8 

723 

3525.50 

4623.58 

81.7 

6.5 

724 

3442.67 

4796.98 

91.7 

6.4 

725 

3478.48 

4818.98 

92.6 

6.4 

726 

3558.35 

4658.77 

83-7 

6.5 

727 

3612.48 

4473.89 

73.4 

6.7 

728 

3635.37 

4279.04 

62.  1 

7.4 

729 

3618.98 

4069. 08 

49.0 

8.8 

730 

3574.84 

3866.91 

35.  1 

11.8 

731 

3587.72 

3809.47 

32.  1 

12.8 

732 

3613.95 

3910.30 

39.  1 

10.7 

733 

3648.48 

4007.29 

46.  1 

9.3 

734 

3676.20 

4121. 41 

53.6 

8.2 

735 

3677.60 

4225.91 

59.8 

7.6 

73  6 

3676.79 

4332.54 

66.0 

7.  1 

737 

3658.87 

4436.08 

71.7 

6.8 

738 

3638. 19 

4545.66 

77.6 

6.6 

739 

3604.38 

4649.50 

83.2 

6.5 

740 

3572.48 

4754.  18 

88.7 

6.4 

74  1 

3514.26 

4836.98 

93.2 

6.4 

78 


RANGE-E 

lAKGE  ACCURACIES  (CONTINUED) 

Fix 

X 

Y 

Angle  of 

Radius  cf 

Ho. 

Coordinate 

Coordinate 

Intersection 

90%  Circle 

742 

3570.50 

4836.50 

92.8 

6.4 

743 

3618.79 

4723. 26 

87.0 

6.4 

744 

3662.90 

4602. 44 

80.9 

6.5 

745 

3691. 10 

4478.20 

74.4 

6.7 

746 

3715.24 

4341.39 

67.  1 

7.0 

747 

3720.  11 

4205. 94 

59.5 

7.5 

748 

3703.21 

4071.  17 

51.4 

8.5 

749 

3678.03 

3939.55 

43.0 

9.8 

75C 

3636. 04 

3801.29 

33.7 

12.3 

751 

3630.00 

3715.02 

28.3 

14.5 

752 

3678.77 

3812. 88 

35.7 

11.6 

753 

3715.73 

3925.56 

43.4 

9.7 

754 

3731.48 

4050.62 

51.0 

8.5 

755 

3756.38 

4168. 71 

58.3 

7.7 

756 

3754.68 

4300. 89 

65.5 

7.  1 

757 

3743.02 

4427.92 

72.2 

6.8 

756 

3710.29 

4554. 17 

78.6 

6.6 

759 

3680.62 

4673.88 

84.6 

6.5 

760 

3627.42 

4792.51 

90.5 

6.4 

76  1 

3572.80 

4899.46 

95.8 

6.5 

762 

3627.  15 

4891.25 

95.1 

6.5 

763 

3683.33 

4778.34 

89.7 

6.4 

764 

3725.77 

4656.01 

83.9 

6.5 

765 

3757. 98 

4525.66 

77.5 

6.6 

766 

3761.56 

4403. 13 

71.5 

6.8 

767 

3798.56 

4272.52 

64.8 

7.2 

768 

3787.95 

4134.55 

57.  1 

7.8 

769 

3772.73 

4001.69 

49.4 

8.7 

770 

3738. 73 

3870.33 

41.  1 

10.2 

771 

3695.62 

3745.45 

32.6 

12.7 

772 

3712.52 

3708.43 

31.2 

13.2 

773 

3777.45 

3829. 70 

40.  2 

10.4 

774 

3803.43 

3957.03 

48.0 

9.0 

775 

3818.28 

4087.24 

55.4 

3.0 

776 

3821.33 

4211. 74 

62.3 

7.4 

777 

3823.70 

4338. 97 

68.8 

6.9 

778 

3813.4  1 

4476. 12 

75.6 

6.7 

779 

3777. 13 

4614.24 

82.  1 

6.5 

780 

3739.58 

4736.39 

87.3 

6.4 

781 

3691.63 

4868. 16 

93.9 

6.5 

782 

3706.65 

4918.98 

96.0 

6.5 

783 

3749.46 

4806. 82 

91.  1 

6.4 

784 

3799.89 

4705. 41 

86.6 

6.4 

785 

3829.30 

4581. 14 

81.0 

6.5 

786 

3857.30 

4455. 00 

75.1 

6.7 

787 

3868.49 

4327.66 

69.0 

7.0 

788 

3872.94 

4198.56 

62.5 

7.3 

789 

3862.59 

4062.95 

55.3 

8.0 

790 

3841. 17 

3928.97 

47.7 

9.0 

791 

3804.62 

3797.  14 

39.4 

10.6 

792 

3759.45 

3596. 86 

27.  1 

15.  1 

793 

3808. 16 

3722.44 

35.  6 

11.6 

794 

3851. 17 

3827.55 

42.7 

9.9 

795 

3879.16 

3942.93 

49.6 

8.7 

796 

3896.63 

4068.89 

56.5 

7.9 

797 

3903.40 

4184. 15 

62.5 

7.3 

798 

3916.60 

4306.  31 

68.8 

6.9 

799 

3904.40 

4428.30 

74.5 

6.7 

800 

3879.66 

4554.  38 

80.2 

6.5 

801 

3851. 19 

4684. 13 

85.9 

6.5 

802 

3810. 16 

4796.95 

90.7 

6.4 

803 

3752.88 

4908. 69 

95.5 

6.5 

804 

3706.66 

4998.96 

99.4 

6.5 

805 

3781. 90 

4945. 89 

97.0 

6.5 

79 


RANGE- 

■RANGE  AC  CUR  AC 

:IES  (CONTINUED) 

Fix 

X 

Y 

Angle  of 

Radius  of 

No. 

Coordinate 

Coordinate 

Intersection 

90%  circle 

806 

3630.54 

4843.53 

92.8 

6.4 

307 

3877.77 

4734.  18 

88.3 

6.4 

808 

3907.61 

4618.  18 

83.4 

6.5 

809 

3932.77 

4495. 95 

78.  1 

6.6 

8  10 

3948. 19 

4376.06 

72.8 

6.8 

811 

3956.79 

4251. 77 

67.0 

7.0 

812 

3953.81 

4125.89 

60.9 

7.5 

813 

3938.32 

3995.79 

54.  1 

8.  1 

814 

3909.94 

3869.78 

46.9 

9.  1 

815 

3864.71 

3744.84 

38.9 

10.7 

816 

3816.98 

3619.07 

30.6 

13.4 

817 

3711.62 

5131.00 

104.5 

6.7 

818 

3788.  16 

5014.61 

99.7 

6.5 

819 

3851.30 

4889.79 

94.7 

6.5 

820 

3901.  13 

4761.08 

89.6 

6.4 

821 

3961.73 

4590.52 

82.7 

6.5 

822 

3977.  10 

4485.23 

78.2 

6.6 

823 

3993.32 

4118.31 

61.5 

7.4 

824 

3996.85 

4193.81 

65.2 

7.2 

825 

3981.72 

4048.49 

57.9 

7.7 

826 

3956.88 

3905.60 

50.2 

8.6 

827 

3916.26 

3758.91 

41.5 

10.  1 

828 

3826.04 

3553.55 

27.6 

14.8 

829 

3961.88 

4691.61 

87.  1 

6.4 

830 

3S07. 91 

4827.37 

92.4 

6.4 

831 

3841.30 

4953.23 

97.3 

6.5 

832 

3789.42 

5085.46 

102.4 

6.6 

833 

3772.12 

5177.05 

105.8 

6.7 

834 

3843.76 

5069. 86 

101.7 

6.6 

835 

3906.65 

4961.46 

97.6 

6.5 

836 

3959.36 

4844.  18 

93.3 

6.4 

837 

4002.25 

4720.39 

88.6 

6.4 

838 

4057. 98 

4659.  81 

86.7 

6.4 

339 

3992.26 

4828.76 

92.8 

6.4   - 

840 

3954. 08 

4915.  16 

96.0 

6.5 

841 

3905.75 

5046.63 

100.8 

6.6 

842 

3822.81 

5167.39 

105.2 

6.7 

843 

3844.43 

5129.36 

103.8 

6.6 

844 

3958.55 

5033.27 

100.3 

6.6 

845 

4016.73 

4919.44 

96.4 

6.5 

846 

4061.27 

4797.13 

92.  1 

6.4 

847 

4159.72 

4597.99 

35.  5 

6.  5 

848 

4119. 01 

4735.07 

90.2 

6.4 

849 

4077.07 

4872.84 

95.0 

6.5 

850 

4021.39 

4990.89 

99.0 

6.5 

851 

3961.68 

5118.03 

103.3 

6.6 

852 

3992.88 

5139.94 

104.0 

6.7 

853 

4056.51 

5026.  15 

100.3 

6.6 

854 

4121. 60 

4874.08 

95.4 

6.5 

855 

4154.57 

4779.44 

92.2 

6.4 

856 

4196.42 

4617.75 

86.7 

6.4 

80 


CLASSIFIED    AZIMUTH-AZIMUTH    POSITIONS 


Contrcl    Stations:       USE    MON    1978   and    MUSSEL    1932 


Standard   Error    Used    in  Computations:    1.3    meters 


Fix 

No. 

619 
620 
621 
622 
623 
624 
625 
626 
627 
628 
631 
632 
633 
634 
635 
636 
637 
638 
639 
640 
641 
642 
643 
644 
645 
646 
647 
648 
649 
650 
651 
652 
653 
654 
655 
656 
657 
658 
659 
660 
661 
662 
663 
664 
665 
666 
667 
668 
669 
670 
671 
672 
673 
674 
675 
677 
678 
679 
680 
681 


Coordinate    Coordinate 


Angle   of 
Intersection 


Radius  of 
90^  circle 


4449.26 

4437.86 

4427.41 

4419.65 

4410.42 

4390.60 

4376.52 

4364.31 

4348.71 

4334.31 

4338.97 

4337.95 

4327.70 

4322.53 

4323. 18 

4324.42 

4327.25 

4394.77 

4386.58 

4377.29 

4367.00 

4355.55 

4345.97 

4256.90 

4260-68 

4264. 77 

4283.59 

4293.24 

4300.30 

4345.53 

4370.38 

4398-77 

4411.48 

4438.30 

4470. 

4502. 

4514. 

4512. 

4520. 

4494. 

4487. 

4477. 

4462. 

4453.74 

4440.56 

4465. 02 

4507.24 

4577.55 

4569.50 

4563.68 

4564.53 

4563.90 

4560.41 

4553.49 

4556.71 

4571.60 

4576.20 

4582.23 

4604.35 

4613.39 


97 
73 
38 
00 
10 
79 
96 
62 
90 


88 
02 
65 
01 


2711.40 

2807.29 

2897.  17 

2989.36 

3072.20 

3154.04 

3234.24 

3319.82 

3411.24 

3502.51 

3381.55 

3290. 68 

3199.92 

3107.83 

3012.49 

2916.00 

2818.39 

2806.86 

2903.09 

2998.63 

3090.26 

3187.06 

3285. 

3516. 

3416. 

3321. 

3208. 82 

3127.36 

3024.70 

3145.58 

3236.09 

3327.77 

3421.46 

3506.23 

3591.21 

3677.50 

3767.26 

3860. 13 

3948.53 

3049.75 

3144. 60 

3243. 86 

3372.74 

3469. 99 

3564. 89 

3652.95 

3743.29 

3595.27 

3681.02 

3776.77 

3872.77 

3965.68 

4057.  10 

4147.06 

4239.64 

4416.22 

4504.92 

4597.03 

4631.33 

4527. 17 


156 

153 

149 

146 

143 

141 

138 

136.2 

133.5 

130.7 

135.  1 

138.7 

142.8 

146. 

151. 

155. 

159 

156 

152 

148. 

145-  1 

141.9 

138. 

133. 

137. 

141. 

144. 

147. 

151. 

144. 

139. 

134. 

130. 

125. 

121. 

117 


5 
6 
4 
1 
8 
7 
8 
1 
1 
2 
2 
9 
6 
4 
114.0 


111 

107. 

139, 

136. 

133. 

129, 

126, 

123. 

119. 

115 

116. 

114. 

111. 

108. 

106, 

103, 

100, 

97, 

92, 

89. 

86 

85.  1 

87.9 


1 

9 
2 
2 
2 
4 
5 
8 
7 
0 
9 
6 
9 
8 
0 
2 
6 
7 
0 
3 
5 


4.  4 
4.2 
4.0 
3.3 

4.  1 
4.4 
4-8 
5.4 

6.  1 

7.  1 
8.7 
7.3 
6.3 
5.6 

5.  1 
4.7 
4.4 
4.0 
4.3 
4.7 
5.  1 
5.5 
6.3 
5.1 
4.5 
4.0 
3.8 
3.5 
3.3 
3.2 
3.  1 
3.0 
2.9 
4.5 
4.2 
4.0 
3.7 
3.5 
3.4 
3.3 
3.  1 
3.2 
3.  1 
3.0 
3.0 
2.9 
2.9 
2.8 
2.8 
2-8 
2.8 
2.8 
2-8 
2.8 
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AZIMUTH- 

■AZIMUTH  ACCURACIES  (CONTINUED] 

i 

Fix 

X 

Y 

Angle  of 

Rad. 

Las  of 

No. 

Coordinate 

Coordinate 

Intersection 

90% 

Circle 

682 

4611.20 

4421.29 

91.0 

2.8 

683 

4610.07 

4315.24 

94.1 

2.8 

684 

4609.84 

4204.66 

97.4 

2.8 

685 

4603. 99 

4098.63 

100.7 

2.8 

686 

4600.52 

3992. 16 

104.0 

2.9 

687 

4601.27 

3883.40 

107.3 

2.9 

688 

4601. 13 

3780.02 

110.  4 

3.0 

689 

4602.22 

3675.49 

113.5 

3.1 

690 

4601.48 

3574.94 

116.5 

3.  1 

691 

4603. 21 

3458.09 

120.0 

3.3 

692 

4524.84 

3728. 78 

114.8 

3.  1 

695 

4657.06 

3506.75 

116.  1 

3.  1 

696 

4648.60 

3602.24 

113.7 

3.  1 

697 

4630. 15 

3696.  13 

111.8 

3.0 

698 

4629.80 

3793.39 

108.9 

3.0 

699 

4623.03 

3889.32 

106.3 

2.9 

700 

4622. 83 

3978.70 

103.7 

2.9 

701 

4617.60 

4071.98 

101.  1 

2.8 

702 

4623.30 

4163.73 

98.2 

2.8 

703 

4618.23 

4256.25 

95.7 

2.8 
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CLASSIFIED    RANGE-AZIMUTH    POSITIONS 
Contrcl    Stations:       MGSSEL    1932    occupied,    initial    USE    MON    1971 
Standard    Errors:      Range — 3    meters;    T-2 — 1.3    meters 


Fix 

X 

Y 

Angle  of 

Radius  of 

No. 

Coordinate 

Coordinate 

Intersection 

90%  Circle 

414 

2898.58 

4611.99 

90.0 

5.  1 

415 

2979.40 

4483.31 

90.0 

5.1 

416 

3027.69 

4342.47 

90.0 

5.  1 

417 

2795.26 

4216.04 

90.0 

5.  1 

418 

3002.90 

4061-  15 

90.0 

5.  1 

419 

2939.58 

3938.67 

90.0 

5.  1 

420 

2904.53 

3888.38 

90.0 

5.  1 

421 

3195.47 

3380.99 

90.0 

5.  1 

422 

3128.22 

3487.82 

90.0 

5.1 

423 

3048. 17 

3599.79 

90.0 

5.  1 

424 

2994.68 

3723.65 

90.0 

5.  1 

425 

2911.61 

3833.20 

90.0 

5.  1 

426 

2833.62 

3921.  16 

90.0 

5.  1 

427 

2749. 17 

4000.71 

9  0.0 

5.  1 

428 

2752.46 

4029.48 

90.0 

5.  1 

429 

2823.93 

3969.92 

90.0 

5.  1 

430 

2894.60 

3920.66 

90.0 

5.  1 

431 

2907.95 

3839.53 

90.0 

5.  1 

432 

2961.44 

3901.64 

90.0 

5.  1 

433 

3005.23 

3980.11 

90.0 

5.  1 

434 

3037.78 

4058.01 

90.0 

5.  1 

435 

3065.09 

4134.  51 

90.0 

5.  1 

436 

3074.88 

4036.74 

90.0 

5.  1 

437 

3C48.02 

3953.93 

90.0 

5.  1 

438 

2996.76 

3871.95 

90.0 

5.  1 

439 

2943. 10 

3809.76 

90.0 

5.  1 

440 

3128.59 

4049. 87 

90.0 

5.  1 

441 

3096.20 

3978.78 

90.0 

5.  1 

442 

3058.84 

3897.24. 

90.0 

5.  1 

443 

3009.61 

3825.86 

90.0 

5.  1 

444 

2914.71 

3837.97 

90.0 

5.  1 

445 

3008.50 

3755.28 

90.0 

5.  1 

446 

3064.52 

3839.64 

90.0 

5.  1 

448 

3153.77 

4023.36 

90.0 

5.  1 

449 

3182.88 

4128.56 

90.0 

5.  1 

450 

3210.38 

4048. 54 

90.0 

5.  1 

451 

3179.45 

3970.85 

90.0 

5.  1 

452 

3146.52 

3889.81 

90.0 

5.  1 

453 

3099.79 

3820.34 

90.0 

5.  1 

454 

3047.75 

3751. 18 

90.0 

5.  1 

456 

3037.99 

3689. 19 

90.0 

5.  1 

457 

3121.90 

3785.62 

90.0 

5.  1 

458 

3186.09 

3883.49 

90.0 

5.  1 

459 

3234.50 

4006.49 

90.0 

5.  1 

460 

3264.92 

4128. 69 

90.0 

5.  1 

46  1 

3079.28 

3672.85 

90.0 

5.  1 

462 

3159.97 

3763.25 

90.0 

5.  1 

463 

3221.69 

3867.47 

90.0 

5.  1 

464 

3273.65 

3994.32 

90.0 

5-  1 

465 

3324.98 

4005.30 

90.0 

5.  1 

466 

3295.53 

3924.46 

90.0 

5.  1 

467 

3257. 07 

3846.60 

90.0 

5.  1 

468 

3218.  14 

3770.87 

90.0 

5.  1 

469 

3172. 68 

3704.85 

90.0 

5.  1 

470 

3122.00 

3646.77 

90.0 

5.  1 

47  1 

3097.64 

3628.62 

90.0 

5.  1 

472 

3089.39 

3573. 64 

90-0 

5.  1 

473 

3118. 69 

3586. 14 

90.0 

5.  1 

474 

3210.37 

3698.21 

90.0 

5.  1 

475 

3262.22 

3765.90 

90.0 

5.  1 
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RANGE-AZ! 

IEUTH    ACCUSACIES     (CONTINUED) 

Fix 

X 

Y 

Angle    of 

Radius   of 

No. 

Coordinate 

Coordinate 

Intersection 

90  *   Circle 

476 

3317.64 

3873. 69 

90.0 

5.  1 

478 

3366.64 

3900.37 

90.0 

5.  1 

479 

3335.51 

3831.40 

90.0 

5.1 

480 

3303. 12 

3760.21 

90-0 

5.  1 

481 

3258.28 

3696. 16 

90.0 

5.  1 

482 

3212.24 

3631.26 

90.0 

5.  1 

483 

3165. 13 

3578.90 

90.0 

5.  1 

484 

3129.45 

3550.77 

90.0 

5.  1 

485 

3142.92 

3499.37 

90.0 

5.  1 

486 

3214.  10 

3575.56 

90.0 

5.  1 

437 

3280.84 

3665.79 

90.0 

5.  1 

488 

3344.31 

3753.54 

90.0 

5.  1 

489 

3395.65 

3847.40 

90.0 

5.  1 

490 

3420.41 

3919.  1  1 

90.0 

5.  1 

491 

3466.94 

3949.20 

90.0 

5.  1 

492 

3439.67 

3856.67 

90.0 

5.  1 

493 

3395.38 

3761.44 

90.0 

5.  1 

494 

3335.71 

3669.35 

90.0 

5.  1 

495 

3275.38 

3583.62 

90.0 

5.  1 

496 

3197.78 

3504. 04 

90.0 

5.  1 

497 

3161.00 

3465. 96 

90.0 

5.  1 
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APPENDIX  C 
PECGEAM  FOB  90  PERCENT  CONFIDENCE  ELLIPSE  PARAMETERS 


C  PRCGF-AM  NAME:  ELLIP 

C 

C  DESCRIPTION:  COMPOTES  ORIENTATION,  90%  SIGMA-X,  90??  SIGMA-Y, 

C  FCR  ERROR  ELLIPSE  ABOUT  A  HYEECGRAPHIC 

C  POSITION  ESTABLISHED  BY  RANGE-RANGE,  AZIMUTH- 

C  AZIMUTH  OR  RANGE-AZIMUTH  POSITION 

C 

C    AUTHCE:     NICHOLAS    E.    PERUGINI 

C  LT.     NOAA 

C  NAVAL  POSTGRADUATE  SCHOOL 

C  DATE:    SEPTEMBER  #  1984 

C 

c 

IMPLICIT  REAL  *  4  (A-H,C-Z) 

C   INITIALIZE  VALUES:   FOR  RANGE-RANGE  AND  AZIMUTH-AZIMUTH: 

C  -XL  AND  YL  ARE  COORDINATES  OF  LEFT  STATION 

C  -XR  AND  YR  ARE  COORDINATES  OF  RIGHT  STATION 

C  -SIGI  AND  SIGR  ARE  RESPECTIVE  STANDAED  EREGES 

C  ASSOCIATED  WITH  EACH  LOP 

C 

C  FOR    RANGE-AZIMUTH: 

C  -XL    ANE    YL    ARE    CCORDINATES    OF    OCCUPIED    STATION 

C  -SIGL    IS    SIGMA    OE    THEODOLITE    LOP 

C  -SIGR    IS    SIGMA    OF    RANGE    LOP 

C******************************************* **************** 

XL=4914.75 
YI=2C09.86 

S3GL=3.0 
C 

XR=2474.75 

YR=4247.42 

SIGR=  3.0 
C 

C  ENTEF  CORRELATION  COEFFICIENT:  USUALLY  ZERO  FOR  R-B, 
C  R-AZ,  AND  AZ-AZ 

RC  =  0.0 

PI=3. 141593 
C 

C  ENTER  INDICATOR  TO  TELL  WHAT  KIND  OF  DATA  IS  ENTERING  PROGRAM 
C  IND  =  1    RANGE-RANGE 

C  IND  =  2    AZIMUTH-AZIMUTH 

C  IND  =  3    RANGE-AZIMUTH 

C***  ************  *******  *******************  *******  ** 

IND=1 
C************** ********  **************************** 

C  INC  IS  A  TOGGLE  WHICH  CHECKS  FOR  BETA  GREATER  THAN  90  DEG. 

C  NOTHING  IN  PROGRAM  SHOULD  BE  CHANGED  FROM  HERE  ON 

C 

C  READ  IN  DATA  FROM  IATA  FILE:  IFIX  =  FIX  NUMBER 

C  FX    =  X  COORDINATE  OF  HYDRO  POSITION 

C  EY    =  Y  COORDINATE  OF  HYDRO  POSITION 

C  TD      ANGLE  OF  INTERSECTION  IN  DEGREES 

C  SENTINEL  IS  IFIX  =  999,  TELLS  PROGRAM  TO  STOP  READING 


C 


10    CONTINUE 
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EEAE  (4,20)  IF2X#ILLePX,FY.TD,RAD 
20    FCRHAT71X#I3,3X,I1,6X,F7.2r6X,F7.2, 5X,F8.4,3X,F5.2) 
INC  =  0 


IF  (IFIX.EC.999)  GO  TO  900 
IP  (ID. LT. 90.)  ED=TD 
IFJID.LT.90.)  GO  TO  30 


C  WORK  WITH  BETA  LESS  THAN  90  DEGREES:  TOGGLE  TURNED  CN  TC  ONE 
DE=180.-TD 
INC=1 
30     CONTINUE 
C 

C  KEEP  TANGENT  FUNCTION  FROM  GCING  UNDEFINED  IN  A  RARE  CASE 
C  OF  TEE  FIX  AND  CONTROL  STATION  HAVING  SAME  COORDINATES 


IF(PX.EQ.XL)  IX=PX+  0.5 
IF(EY.EQ.YL)  PX=PY+  0.5 


C 

C  CHANGE  EEGREES  TO  EADIANS 

3ETA=.01745329*DD 
C  USE  LEFT  STATION  AS  BASIS  FCR  COMPUTATIONS 

C  ORIENTATION  ANGLES  WILL  BE  FIXED  WITH  RESPECT  TO  LEFT  LOP 
C 

c 

C  FINE  AZIMUTH  FROM  NCRTH  BETWEEN  HYDRO  POSITION  AND  LEFT 
C  STATION.  AZIMUTH  WILL  EE  DEFINED  BETWEEN  0-180  EEGREES 
C  MEASURED  CLOCKWISE  EROM  NORTH. 

C  THIS  IS  THE  RANGE-RANGE  AZIMUTH  DETERMINATION. 
IF  (INE.NE.  1)  GO  TO  40 
IF(PY.GE.YL)  TEEN 

IF(PX.GE.XL)     THEN 

ALPHA    =    PI-ATAN  (  (PY-YL) /  (PX-XL) ) 
ELSE 

ALPHA   =       ATAN  (  (PY-YL) /( XL-PX)  ) 


EISE 


C 


END    IF 

IF(PX.GE.XL)     THEN 

ALPHA   =    ATAN  (  (YL-PY)  /  (PX- XL)  ) 

ELSE 

ALPHA   =    PI-ATAN ( ( YL-PY) / ( XL-PX)  ) 

END    IF 
END    IF 
GC    TO    60 
:    AZIMUTH    FIXING    FOR    AZIMUTH-AZIMUTH    POSITIONS 
40  CONTINUE 

IF(EY-GE.YL)     TEEN 

IF(PX.GE.XL)     THEN 

ALPHA   =    ATAN     ( (PX-XL) / (PY- YL) ) 
ELS  E 

ALPHA   =    PI-ATAN  (  (XL-PX) /(PY-YL)) 
END    IF 
EISE 

IF(PX.GE.XL)     THEN 

ALPHA  =    PI-ATAN  (  (PX-XL)/ (YL-PY)  ) 
ELSE 

ALPHA    =    ATAN ( (XL-PX) /(YL-PY) ) 
END    IF 
END    IF 
C 

C       AZIMUTH    EQUALS    THEIA    FOR    RANGE    AZIMUTH    CASE,     ASSUMING 
C       THEODOLITE    SIGMA    IS    LESS    THAN    RANGE    SIGMA 
C 

3F     (IND.EQ.3)     GC    TO    70 

C 

C    BEGIN    COMPUTING    THEIA,     THAT    IS    THE    ANGLE    OF    ROTATION     FROM 
C    LEFT    LCP 
60  CONTINUE 

B1=SIGL**2*SIN  (2*BETA)  +2*RO*SIGL*SIGR*SIN ( BETA) 
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E2=SIGL**2*COS  (2*BETA)  +2*R0*SIGL*SIGR*C0S  (BETA)  4SIGR**2 

IE  (ABSJB2) .LI. 0.000 1   E2=.0001 

E2=E1/B2 
C  COMEUTE  ROTATION  ANGLE  FROM  LEFT  LOP 

IH=0.5*ATAN (B3) 
90   CCNTINOE 
C 

c 

C    DEEINE    SEMI-MAJOR    AXES    ORIENTATION    IN    TERMS    OF    0-180    EEGREES 

C    ROTATION,     CLOCKWISE    FROM    NOETH 

C 

C    RANGE-RANGE    CASE 

IF(IND.EQ.I)  TEEN 

IF(INC-EQ.I)  THETA=AIPHA+TH 
IFjINC.EQ.O)  THETA=AIPHA-TH 
END  IF 
C 
C  AZIMUTE-AZIMUTH  CASE 

IF  (IND.EQ.2)  TEEN 

IF{INC.EQ.C)  THETA=AIPHA+TH 
IF(INC.EQ.I)  THETA=AIPflA-TH 
END  IF 
C 

c 

C    FIX    ROTATION    ANGLE    EROM    0-180    DGREES 

70  CCNTINOE 

C    CONDITION    FOR    RANGE- AZIMOT H    DATA 
II(IND.EQ.3)     TEETA=ALPHA 
IF  (THETA.LT.O.)     THETA=    EI+THETA 
IEJTHETA.GT.PI)     THETA=    THETA-PI 
C    DEG    IS    THE    SEMI-MAJCE    ELLIPSE    AXIS    ORIENTATION     IN    DEGREES 

DEG=57. 295779*1HETA 
C    COMEUTE    90%    SIGMAX    AND    SIGMAY    OF    ERROR    ELLIPSE 

CALL    PROB  (SIGI.SIGB,BO.TD,SGX90fSGY90,RADIUS#ELAE#CIRAR) 
WEITE(7-  100)        IFIXjPX.PY.TD-SGXSO.SGYSiOfDEG 
100       FCEMATflX,l3,3X,F7.2,3X,F7.2,3X,F5.  1 ,  3X,  F  4.  1 ,  3X  ,  F4  . 1  , 
*    3X,F5.1) 
GC    TO    10 
900      CCNTINOE 
STOE 
END 
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APPENDIX    D 
ACCURACY    CLASSIFICATION:    90    PERCENT    CONFIDENCE    ELLIPSIS 

CLASSIFIED    RANGE-RANGE    POSITIONS 
Contrcl    Stations:      BI2CH   LAB    1982    and   MUSSEL    1932 
Standard    Error    Used    ir   Computations:    3    meters 


Fix 

X 

Y 

Beta 

90% 

30% 

Orienta 

No. 

Coord. 

Coord. 

Si^ma   X 

Sigma    Y 

tion 

1 

2668.05 

4942.07 

127.0 

10.2 

5.  1 

79.0 

2 

2852.01 

5076.69 

121.6 

9.3 

5.2 

85.3 

3 

3041.58 

5194.85 

1  18.6 

8.9 

5.3 

90.2 

4 

3040.27 

5148.74 

117.0 

8.7 

5.3 

90.6 

5 

2838.26 

5021.35 

120.3 

9.1 

5.2 

85.3 

6 

2640.06 

4877.13 

126.9 

10.2 

5.1 

78.1 

7 

2553.46 

4752.36 

130.4 

10.9 

5.0 

74.1 

8 

2724.10 

4885.47 

121.4 

9.3 

5.2 

82.0 

9 

2889.06 

5013.63 

1  17.  6 

8.8 

5.3 

37.2 

10 

3C75.61 

5124.37 

1  15.0 

8.5 

5.4 

91.9 

11 

3172.63 

5136.23 

1  12.7 

8.2 

5.5 

94.5 

12 

2958.  13 

5006.65 

1  14.4 

8.4 

5.4 

89.7 

13 

2771.33 

4876.02 

1  18.0 

8.8 

5.3 

84.2 

14 

2581.02 

4729.89 

127.0 

10.2 

5.1 

75.9 

15 

2584.41 

4665.31 

124.0 

9.7 

5.2 

76.7 

16 

2140.05 

4805.63 

1  16.7 

8.7 

5.3 

33.8 

17 

2913.27 

4928.51 

1  12.8 

8.2 

5.5 

89.2 

18 

3097.44 

5047.77 

111.2 

8.1 

5.5 

93.5 

19 

3193.94 

5103.46 

1  10.9 

8.0 

5.5 

95.5 

20 

2627.47 

4659-05 

1  18.8 

8.9 

5.3 

79.8 

21 

2762.15 

4769.60 

113.2 

8.3 

5.5 

85.4 

22 

2904.91 

4875.24 

1  10.5 

8.0 

5.5 

89.7 

23 

3056.38 

4975.73 

109.3 

7.9 

5.6 

93.3 

24 

3207.36 

5064.00 

108.9 

7.8 

5.6 

96.3 

25 

3373.36 

5101.96 

107.  1 

7.7 

5.7 

100.0 

26 

2190.64 

5007.38 

106.8 

7.6 

5.7 

96.7 

27 

3015.45 

4900.75 

107.1 

7.7 

5.7 

93.2 

28 

2839.38 

4778.31 

108.7 

7.8 

5.6 

88.3 

29 

2679.99 

4651.31 

1  12.8 

8.2 

5.5 

83.4 

30 

2727.88 

4632.77 

106.9 

7.6 

5.7 

86.7 

31 

2899.83 

4776.97 

105.2 

7.5 

5.7 

91.3 

32 

3092.77 

4901.75 

104.4 

7.4 

5.8 

95.6 

33 

3295.42 

5019.92 

105.0 

7.5 

5.7 

99.2 

34 

3502.07 

5123.36 

106.0 

7.6 

5.7 

102.6 

35 

3697.02 

5151.09 

105.3 

7.5 

5.7 

106.2 

36 

3474.77 

5065.57 

104.  1 

7.4 

5.8 

102.7 

37 

3257.48 

4953.37 

102.7 

7.3 

5.8 

99.3 

38 

3043.03 

4821.88 

101.7 

7.2 

5.9 

95.5 

39 

2845.07 

4680.80 

101.7 

7.2 

5.9 

91.4 

40 

2746.96 

4608.32 

103.  1 

7.3 

5.8 

88.6 
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CLASSIFIEE    AZIMUTH-AZIMUTH    POSITIONS 
Control    Stations:       USE   MON    1978   and    MUSSEL    1932 
Standard   Error    Used    in   Computations:    1.3    meters 


Fix 

X 

Y 

Beta 

90% 

90fl 

Orienta- 

No. 

Ccord. 

Coord. 

Sigma  X 

Sigma  Y 

tion 

619 

4449.26 

271  1.40 

156.9 

9.8 

2.0 

139.  4 

620 

4437.86 

2807.29 

153.0 

8.5 

2.0 

139. S 

621 

iJ £427-  U 1 

2897.17 

149.6 

7.5 

2.0 

139.8 

622 

4419.65 

2989.36 

146.2 

6.8 

2.  1 

139.8 

623 

441C.42 

3072.20 

143.4 

6.3 

2.1 

139-6 

624 

4390.60 

3154.04 

141.3 

5.9 

2.1 

139.  1 

625 

4376.52 

3234.24 

138.9 

5.6 

2.  1 

138.6 

626 

4364.31 

3319.82 

136.2 

5.3 

2.  1 

138. C 

627 

4348.71 

341  1.24 

133.5 

5.0 

2.1 

137.3 

628 

4334.31 

3502.51 

130.7 

4.7 

2.2 

136.5 

631 

4338.97 

3381.55 

135.1 

5.2 

2.1 

137.4 

63  2 

4337.95 

3290.68 

138.7 

5.6 

2.1 

137.8 

633 

4327.70 

3199.92 

142.8 

6.2 

2.1 

138.  1 

634 

4322.53 

3107.83 

146.9 

6.9 

2.  1 

138.2 

635 

4323-18 

3012.49 

151.0 

7.9 

2.0 

138.3 

636 

4324.42 

2916.00 

155.3 

9.2 

2.0 

138.  1 

637 

4327.25 

2818.39 

159.8 

11.3 

2.0 

137.7 

638 

4394.77 

2806.86 

156.0 

9.5 

2.0 

138.9 

639 

4386.58 

2903.09 

152.0 

8.2 

2.0 

139.  1 

640 

4377.29 

2998.63 

148.4 

7.2 

2.  1 

139.  1 

641 

4367.00 

309C.26 

145-1 

6.6 

2.1 

138.9 

642 

4355.55 

3187.06 

141.9 

6.0 

2.  1 

133.5 

643 

4345.97 

3285.88 

138.5 

5.6 

2.1 

138.0 

644 

4256.90 

3516.02 

133.6 

5.0 

2.1 

135.5 

645 

4260.68 

3416.65 

137.4 

5.4 

2.1 

136.  2 

646 

4264.77 

332  1.01 

141.  1 

5.9 

2.1 

136.8 

647 

4283.59 

3208.82 

144.8 

6.5 

2.1 

137.5 

648 

4293.24 

3127.36 

147.7 

7.  1 

2.  1 

137.8 

649 

4300.30 

3024.70 

151.8 

8.1 

2.0 

137.9 

650 

4345.53 

3145.58 

144.  1 

6.4 

2.  1 

138.5 

651 

4370.38 

3236.09 

139.1 

5.7 

2.1 

138.5 

652 

4398.77 

3327.77 

134.2 

5.  1 

2.1 

138.  4 

653 

441  1.48 

3421.46 

130.2 

4.7 

2.2 

138.0 

654 

4438.30 

3506.23 

125.9 

4.3 

2.2 

137.7 

655 

4470.97 

3591.21 

121.6 

4.0 

2.3 

137.4 

656 

4502.73 

3677.50 

117.4 

3.8 

2.3 

137. C 

657 

4514.38 

3767.26 

1  14.0 

3.6 

2.4 

136.2 

658 

4512.00 

3860.13 

111.1 

3.5 

2.4 

135.2 

659 

4520.  10 

3948.53 

107.9 

3.4 

2.4 

134.3 

660 

44S4.79 

3049.75 

139.2 

5.7 

2.  1 

141. C 
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CLASSIFIED    RANGE-AZIMUTH    POSITIONS 


Control    Stations 

:       USE    MON 

1978, 

MUSSEL     ' 

I932 

Stand 

ard   Errors: 

RANGE — 3 

meters; 

T-2—  1. 

,3 

Fix 

X 

Y 

3eta 

9  0% 

90% 

Orierta 

No. 

Coord. 

Coord. 

Sigma    X 

Sigma   Y 

tion 

414 

2898.58 

461  1.99 

90.0 

6.4 

2.8 

49.3 

415 

2979.40 

4483.31 

90.0 

6.4 

2.8 

64.9 

416 

3027.69 

4342.47 

90.0 

6.4 

2.8 

80.2 

417 

2795.26 

4216.04 

90.0 

6.4 

2.8 

95.6 

418 

3CC2.90 

4061.15 

90.0 

6.4 

2.8 

1C9.4 

419 

2939.58 

3938.67 

90.0 

6.4 

2.8 

123.6 

420 

2904.53 

3888.38 

90.0 

6.4 

2.8 

129.9 

421 

3  195.47 

3380.99 

90.0 

6.4 

2.8 

140.2 

422 

3128.22 

3487.82 

90.0 

6.4 

2.8 

139.3 

423 

3048.  17 

3599.79 

90.0 

6.4 

2.8 

138.5 

424 

2994.68 

3723.65 

90.0 

6.4 

2.8 

135.2 

425 

2911.61 

3833.20 

90.0 

6.4 

2.8 

133.5 

426 

2833.62 

3921.16 

90.0 

6.  4 

2.8 

132.3 

427 

2749.  17 

400C.71 

90.0 

6.4 

2.8 

132.0 

428 

2752.46 

4029.48 

90.0 

6.4 

2.8 

128.  1 

429 

2823.93 

3969.92 

90.0 

6.4 

2.8 

128.  5 

430 

2894.60 

3920.66 

90.0 

6.4 

2.8 

127.9 

431 

2907.95 

3839.53 

90.0 

6.4 

2.8 

133.3 

432 

2961.44 

3901.64 

90.0 

6.  4 

2.8 

125.4 

433 

3C05.23 

3980.11 

90.0 

6.4 

2.8 

116.7 

434 

3C37.78 

4058.01 

90.0 

6.4 

2.8 

103.6 

435 

3C65.09 

4134.51 

90.0 

6.4 

2.8 

100.  8 

436 

3C74.88 

4036.74 

90.0 

6.4 

2.8 

1C9.3 

437 

3C48.02 

3953.93 

90.0 

6.4 

2.8 

117.  1 

438 

2996.76 

3871.95 

90.0 

6.4 

2.8 

125.7 

439 

2943. 10 

3809.76 

90.0 

6.4 

2.8 

133.  1 

440 

3128.59 

4049.87 

90.0 

6.4 

2.8 

106.  8 

44  1 

3096.20 

3976.78 

90.0 

6.4 

2.8 

113.  4 

442 

3C58.84 

3897.24 

90.0 

6.4 

2.8 

120.9 

443 

3009.61 

3825.86 

90.0 

6.4 

2.8 

128.2 

444 

2914.71 

3837.97 

90.0 

6.4 

2.8 

132.9 

445 

30C8.50 

3755.28 

90.0 

6.4 

2.8 

132.7 

446 

3C64.52 

3839.64 

90.0 

6.4 

2.8 

124.7 

448 

3153.77 

4023.36 

90.0 

6.4 

2.8 

108.3 

449 

3182.88 

4128.56 

90.0 

6.4 

2.8 

99.5 

450 

3210.38 

4048.54 

90.0 

6.4 

2.8 

105.  1 

451 

3179.45 

3970.85 

90.0 

6.4 

2.8 

111.4 

452 

3146.52 

3889.81 

90.0 

6.4 

2.8 

1  18.0 

453 

3C99.79 

3820.34 

90.0 

6.4 

2.8 

124.3 

454 

3047.75 

3751.18 

90.0 

6.4 

2.8 

130.9 

90 


IIST  OF  REFERENCES 

1.  Kaplan,  A.,  Error  Analysis  of  Hydroqraphic  P  ositioning 
and  the  Application  of  LeasT  Squares'^  M.  S.  TKesis, 
flaval  "Postgra^uaTe  School,- Honterey,  1980. 

2.  Heinzen,  M.  R.<  Hydrograrhic  Surveys:   Geodetic 
Ccntrcl  Criteria,  Ml  5.  Thesis, "CornellHTnTver sit y, 

3.  Naval  Warfare  Research  Center  Research  Memorandum 
NWRC-RM34,  Stanford  Research  Institute,  Menlo  Park, 
CA.,  Mathematical  Considerations  Pertaining  to  the 
Accuracy  of  Position  ■Location  and  "NaviqaTTTon  "SysTems, 
l$.rZ   TToy  W.  "3.  Bart  ana"  oTEers,  T9"5T>. 

4.  Adams,  K.  T. ,  Hydroqraphic  Manual,  U.S.  Coast  and 
Geodetic  Survey,  1942. 

5.  International  Hjdrographic  Organization  Special 
Pullication  No.  44  (2nd  Edition),  International 
Hydrographic  Organization  Standards  for  "Hyclro graphic 
Surveys  anq*  ClassiTicaTicn  grifefia  "for  Deep  Sea 
Soundings,  19H'. 

6.  Wallace,  J.  L.  ,  telephone  conversation,  National 
Oceanic  and  Atmospheric  Administration,  National  Ccean 
Service,  Rockville,  MD,  June  1984. 

7.  Laurila,  S.  H.  ,  Electronic  Surveying  and  Navigation , 
John  Wiley  &  Sens,  T97"5. 

8.  Davis.  R.  E.  ,  and  others,  Surveying  Theory  and 
Practice,  McGraw-Hill,  1981. 

9.  Munson,  R.  C. ,  Positioning  Systems:   Report  on  the 
Work  of  WG414b,  paper  presented  at  the  XT- 

Infer national  congress  of  Surveyors,  Stockholm, 
Sweden,  1977. 

10.  Bowditch,  N. ,  American  Practical  Navigator,  Defense 
Mapping  Agency- EycTE ographTc  Cenler,  19  77. 

11.  Hoover,  W.  E.,  Algorithms  for  Confidence  Circles  and 
Ellipses,  unpublished  paper,  National  Oceanic  and 
Atmospheric  Administration,  Woods  Hole,  MA,  June  1984. 

12.  ACIC  Technical  Beport  No.  96,  USAF  Aeronautical  Chart 
and  Information  Center,  St.  Louis.  MO.  Principles  of 
Error  Theory  and  Cartoqraphic  Applications,  fy  C.  R". 
Green waTE~and  AT  E.  Schultz,  19627 


91 


13.  Waltz,    D.    A.,    An    Error    Analysis   of    Range- Azimuth 

Positioning,    M.    S.^TEesis ,    Naval    Postgraduate    School, 
flcnTerey,    1983. 


14. 


Umbach,    M.    J. ,    Eydroqraphic    Manual    Fourth   Edition 
National   Oceanic   and   Atmospheric    Administration, 
National   Ocean    Survey,    1976. 


15.  Harter,    H.    L. ,    "Circular    Error    Probabilities,"    Journal 
of   the   American  Statistical   Association,    Vol.    53, 
December    T96TT. 

16.  NPS    Technical    Ncte    VM-.12,    Osing   DISSPLA   at    NPS,    Naval 
Postgraduate    School,    Monterey,    June   T9~54. 

17.  National   Ocean    Survey    Technical   Manual   No.     1., 
HYEEOPLOT^HYDRCIOG    SYSTEM    MANUAL,     by    J.    L.     Wallace, 


92 


INITIAL  DISTRIBUTION  LIST 


No, 


3. 


4. 


6. 


7. 


10, 


Defense  Tec 
Cameron  Sta 
Alexandria, 


lihxar' 
Naval  1 
Monterey 


Co 

fa 


hnical  Information  Center 
tion 
VA   22314 

de  0  1^42 

raduate  School 
A   92943 


Copies 
2 


Chairman,  Department  of  Oceanography 
Code  68Mr 


Naval  Post 
Monterey 


fa 


NCAA  liaiso 
Post  Office 
Monterey,  C 


raduate  School 
A   93943 

n  Officer 
Box  8688 
A   93943 


Erofessor  Joseph  J.  von  Schwind 
Code  68Ys 
Naval  Post< 
Monterey, 


raduate  School 
A   93943 


It.  Nichola 
NOAA  Ship  E 
1801  Fairvi 
Seattle,  WA 

Cffice  of  t 

Naval  Ocean 

Department 

Washington, 

Ccmmander 
Naval  Ocean 
NSTI,  MS   3 

Ccmmanding 
Naval  Ocean 
Eay  St.  Lou 
NSTI,  MS   3 

Commanding 
Naval  Ocean 
Eay  St.  Lou 

NSTI,    MS      3 


s   E.    Eerugini 
AINIEE 

ew   Avenue,    East 
98  1C2 

he  Director 

ography  Division  (OP-952) 
of  the  Navy 
D.C.   20350 


ography  Command 
9529 

Officer 

ografhic  Office 
is 
9522 

Officer 

Research  and  Development  Activity 
is 
9522 


93 


11.   Chief  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA   21217 


12.  Chairman,  Oceanccraphy  Department 
D.  £.  Naval  Academy 
Annapolis,  MD   2  1402 

13.  Chief,  Hydrographic  Programs  Division 
Defense  flapping  Agency  (Code  PPH) 
Euilding  56 

C.S.  Naval  Observatory 
Washington,  D.C.   20305-30C0 

14.  Associate  Deputy  Director  for  Hydrography 
Defense  Mapping  Agency  (Cede  DH) 
Euilding  56 

C.S.  Naval  Observatory 
Washington,  D.C.   20305-3000 

15.  Topographic  Sciences  Department 
Code  TSD-MC 

Defense  Happing  School 

Ft.  Belvoir,  VA   22060-5828 

16.  Director  Charting  and  Geodetic  Services 
N/CG,  Room  1006,  WSC-1 

National  Oceanic  and  Atmospheric  Administration 
Eockville,  MD   2C852 

17.  Chief,  Nautical  Charting  Division 
N/CG2,  Room  1003,  WSC-1 

National  Oceanic  and  Atmospheric  Adminsitrat ion 
Eockville,  MD   20852 

18.  Chief,  Kydrographic  Surveys  Branch 
N/CG24,  Room  404,  WSC-1 

National  Oceanic  and  Atmospheric  Administration 
Eockville,  MD   20852 

19.  Program  Planning,  Liaison,  and  Training  Div. 
NC2,  Room  105,  Rockwall  Building 

National  Oceanic  and  Atmospheric  Administration 
Eockville,  MD   20852 

20.  Director,  Pacific  Marine  Center 
N/MCP 

National  Ocean  Service,  NCAA 
1801  Fairview  Avenue,  East 
Seattle,  WA   98  102 

21.  Director,  Atlantic  Marine  Center 
N/MCA 

National  Ocean  Service,  NCAA 
439  West  York  Street 
Norfolk,  VA   23510 

22.  Chief.  Program  Services  Division 
N/MCA2 

National  Ocean  Service,  NOAA 
439  West  York  Street 
Norfolk,  VA   23510 

23.  It.  Bruce  Hillard.  NOAA 
NOAA  Ship  DAVIDSCN 

1801  Fairview  Avenue,  East 
Seattle,  WA   98 1C2 


94 


24.  Mr.  Mark  Faye 
430  Cedar  Lane 

last  Meadow,  NY   11554 

25.  Mr.  John  Kolesar 
278  Tripp  ST. 
Swoyersville,  P3   18704 

26.  Mr.  Edward  N.  Perugini 
72  Maltby  Ave. 
Swoyersville,  PS   18704 


95 


13      3  7 


210521 


Perugini 

Criteria  for  the 
classification  of 
hydrographic  position-* 
ing  data. 


.'4   JUN    66 


3  3218 


21C521 


Thesis 

P363 

c.l 


Perugini 

Criteria   for  the 
classification  of 
hydrographic  position- 
ing data. 


